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Dietz,  James  Eric 
N.S.  in  Chemical  Engineering 
Rose-Hulman  Institute  of  Technology 
January  1986 

Economic  Comparison  of  a  High  Pressure  Crystallization 
Process  with  a  Distillation  Process 
Dr.  W.  W.  Bowden 

A  process  capable  of  a  99-*-%  purity  was  desired. 
There  exists  distillation  technology  capable  of  this 
separation.  But,  the  heat  added  and  removed  in  a 
distillation  process  is  very  high  for  a  high  purity 
operation  and  is  very  expensive.  A  more  energy  efficient 
process  is  desired  by  using  a  high  pressure  crystallization 
process.  Crystallization  of  many  systems  may  result  in  a 
complete  separation  and  at  theoretical  yields.  There  were 
two  designs  examined  for  high  pressure  crystallization.  Th 
designs  differed  in  the  pressurization  of  a  batch  charge. 
The  first  was  pressurized  by  an  inert  gas.  The  gas  was 
compressed  and  fed  to  the  charged  crystallizer.  The  second 
was  hydraulically  pressurized.  This  crystallizer  was  a 
double  sided  piston  and  cylinder  arrangement.  The  low 
pressure  end  had  a  large  piston  area  and  the  high  pressure 
end  had  a  smaller  area.  The  pressure  difference  being  the 


ratio  of  the  areas.  It  was  found  that  the  hydraulic  was 
least  costly,  $  171,300,  to  build  and  operate  on  a  10  year 
basis.  The  gas  pressurized  process  was  next  at  $  180,795. 
The  distillation  process,  as  expected,  was  most  costly  at 
$  196,200.  The  gas  pressurized  process  may  have  more 
promise  to  further  reduce  operating  costs.  The  inert  gas 
could  be  replaced  with  air  and  a  diaphragm  could  be  used  to 
separate  it  for  the  liquid  charge.  This  would  result  in  an 
annual  savings  of  $  4,540  and  a  savings  of  $  6,600  in 
capital  expense.  The  gas  pressurization  process  is  also  the 
process  on  which  a  patent  would  most  probably  be  obtained. 
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I.  Introduction 


Crystallization  is  the  formation  of  solids  from 

solutions  or  melts.  This  process  may  be  accomplished  in 

three  methods.  The  first  involves  separation  through  a 

vapor-solid  equilibrium.  Generally,  this  involves  a 

solid  which  is  slightly  soluble  in  the  gas  phase,  for 

example,  iodine  crystals.  The  crystals  sublime  and 

reform  on  some  new  site.  This  has  been  called  fractional 

sublimation.  Second  is  a  solid-liquid  equilibrium 

process.  Commonly,  a  solid  can  be  purified  by  dissolving 

it  in  a  liquid  with  which  it  is  slightly  soluble  at  room 

temperature  and  very  soluble  at  some  higher  temperature. 

It  is  assumed  that  the  impurity  present  in  small 

quantities  is  not  favored  to  deposit  in  the 

recystallization.  This  is  called  recrystallization  or 

solvent  recrystallization.  The  third  method  is  called 

fractional  solidification  or  melt  crystallization.  This 

method  works  best  and  is  usually  used  with  concentrated 

systems  with  an  objective  of  very  high  purity.  The  last 

method  is  of  interest  due  high  recovery  and  absence  of 

solvent,  another  contaminant.  In  the  melt,  the  substance 

itself  acts  as  the  solvent,  eliminating  the  need  of  a 
2  fl 

solvent.  Melt  crystallization  usually  has  the 
advantage  of  lower  operating  costs. 

Melt  crystallization  is  not  a  recent  discovery.  It 


2 


was  first  used  during  the  nineteenth  century  by  whalers. 

They  found  that  potable  water  could  be  obtained  from  sea 

ice  or  icebergs  if  melted  properly.^®  Also  by  1918  crude 

paraffin  waxes  were  separated  from  petroleum  oils 

exclusively  by  crystallization.  These  waxes  were  further 

refined  by  freezing  the  wax  in  sheets  and  placing  them 

over  heating  coils.  Slow  melting  proceeded,  resulting 

first  in  the  appearance  of  oil  beads  which  were  allowed 

29 

to  run  off  leaving  an  even  more  pure  product.  The 

first  industrial  crystallizer,  the  Proabd  refiner,  was 

developed  in  1960.  This  batch  process  slowly  cooled  a 

charge  until  solid  and  then  slowly  reheated  it.  The 

impurities  melt  and  drain  first  and  leave  a  progressively 

20 

more  pure  product.  This  work  will  proceed  further  than 
the  Proabd  refiner.  It  uses  high  pressure,  rather  than 
temperature,  as  a  driving  force.  The  high  pressure 
process  will  be  less  costly,  while  achieving  higher 
purity  than  a  distillation  process  in  use  today. 


The  stream  to  be  refined  has  a  composition  of  90% 
benzene  and  a  desired  product  composition  of  99-»-% 
benzene . 

1.  Design  two  high  pressure  crystallization 
processes,  gas  pressurization  and  hydraulic 
pressurization,  capable  of  the  desired  separation. 

2.  Design  a  distillation  process  capable  of  the 
desired  separation. 

3.  Compare  the  three  designs  on  capital  investment 
and  on  utility  costs. 

4.  Write  a  Flowtran  compatible  program  to  design  the 
crystallizer  with  gas  pressurization. 
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III .  Literature  Summary 


Crystallization 


4 


For  discussion  of  crystallization,  there  must  be 
several  preliminary  topics  discussed.  First  nucleation 
must  be  discussed  to  understand  the  formation  of 
crystals.  Next  phase  diagrams  of  solid-liquid  systems 
must  be  understood.  Also,  for  design  of  a  new  process, 
the  existing  processes  must  be  reviewed.  Finally,  the 
energy  requirement  of  each  processes  must  be  determined. 


....... 
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whatever  the  driving  force  present,  nucleation 
proceeds  in  three  fundamental  steps.  The  first  step 
demands  that  the  solution  obtain  some  degree  of 
supersaturation  or  subcooling.  Supersaturation  is 
cooling  beyond  the  normal  melting  point  of  the  system. 
Supersaturation  is  not  sufficient  cause  for  the 
development  of  crystals.  Microscopic  seeds  or  nuclei  must 
be  present  for  each  successive  molecule  to  deposit. 

These  seeds  may  be  small  crystals  or  any  kind  of 


microscopic  debris  or  dust.  The  third  step  is  the  growth 

of  the  crystal  around  the  developed  centers.  Figure  1 

will  help  explain  some  other  terms  describing 

28 

crystallization. 


NUCLEATION 


SECONDARY 
(crystal  induced  * 


HOMOGENEOUS  HETEROGENEOUS 
( spontaneous ) 

Figure  1.  Hierarchy  of  Crystallization  Terms 

Most  crystals  are  formed  around  some  sort  of  dust  or 

debris,  heterogeneous  nucleation.  Most  reported  cases  of 

spontaneous  crystallization  are  in  fact  seeded  in  some 

manner  such  as  by  atmosphereic  dust.  Normal  aqueous 

6  8 

solutions  contain  10  -  10  particles  per  cubic 

2  8 

centimeter.  Industrially,  crystallization  is  carried 
out  through  heterogeneous  or  secondary  nucleation. 
Deliberate  seeding  gives  some  control  over  product  size 
and  distribution.  The  rainmaking  silver  iodide  is  an 
example.  It  has  the  ability  to  crystallize  water  in  the 
atmosphere  due  to  the  similarity  of  the  Agl  crystal 
lattice  and  that  of  H2O.  Many  factors  are  involved  in 
the  ability  of  a  material  to  seed  crystals.  Lattice 


PRIMARY 
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similarity  may  or  may  not  contribute  to  the 

26 

crystallization  process. 


Phase  Diagrams 


Knowledge  of  phase  diagrams  is  required  to 
understand  the  operation  of  a  crystallizer.  The  phase 
diagram  is  the  equilibrium  relationship  at  specified 
concentration  and  one  other  extensive  variable  (pressure 
or  temperature).  This  work  focuses  on  the  pressure 
variation . 

For  discussion,  an  example  of  a  binary  system  with 
eutectic  and  continuous  solid  solution  will  be  examined 
in  Figure  2.  and  P2  are  pressures  corresponding  to 

pure  component  transition  pressures.  At  these  points  the 
liquid,  mother  liquor,  and  solid  of  these  pure  components 
are  at  equilibrium.  The  point,  e,  at  the  intersection  of 
the  two  curves  is  called  the  eutectic  point.  The  two 
curves  go  through  e  due  to  the  effect  of  subcooling.  The 
pressure  increase  or  supersaturation  is  required  to  form 
crystals.  The  pressure  corresponding  to  this  point  is 


I  i 


Mi 
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Now  that  the  basic  points  of  interest  on  the  phase 


diagram  have  been  examined,  further  discussion  may 
proceed  on  the  usefulness  of  phase  diagrams.  Assuming  a 
solution  of  composition  X^,  the  pressure  can  be  increased 
to  where  the  first  solid  component  2  should  begin  to 
form.  As  the  amount  of  solid  grows,  the  concentration  of 
component  1  grows  in  the  mother  liquor.  At  any  desired 
separation,  there  are  two  phase  concentrations.^^ 

Figure  2  can  be  used  to  find  the  relative 
quantities  of  solid  and  liquid  using  mass  balances.  This 
known  as  the  lever  rule,  involves  measuring  certain 
lengths  on  the  graph. 


Mass  balances 

Overall:  F  *  1  +  s 

Component  2:  *  ^'^1  *  ®'^s 

F-X^  =■  1-X^  +  (F-1)  X^ 


distance  between  X„  and  X„ 

1/F  .  £ - § 

distance  between  X,  and  X 
Also:  ^  ^ 

distance  between  X,  and  X 

s/r  - - 1 - £- 

distance  between  X,  and  X 

i  s 


(1) 

(2) 

(3) 

(4) 

(5) 


(6) 


Binary  mixtures  to  some  degree  do  not  crystallize 


out  in  pure  form.  Due  to  an  inconsistent  crystal 


lattice,  some  impurity  may  be  introduced.  The  result  can 
be  a  solid  solution.  Some  are  completely  miscible  but 
others  may  have  very  small  solubilities. 

There  are  other  relationships  in  addition  those 
shown  on  Figure  2  which  also  must  be  understood  to 
discuss  high  pressure  crystallization.  The  first  of 
these  is  the  relationship  between  pressure  and  volume 
decrease  under  isothermal  conditions.  Figure  3. 

Curve  A  represents  the  phase  change  for  the  pure 
liquid  material.  Crystallization  proceeds  by  compression 
at  constant  pressure,  to  A2,  beyond  A, 
all  material  is  solid  and  compression  of  the  solid  phase 
again  marks  a  rise  in  pressure.  The  curve  B  represents 
the  impure  material.  The  difference  between  Aj^  and  is 
the  transition  pressure  elevation.  Like  the  freezing 
temperature  depression,  the  pressure  difference  is  a 
function  of  the  amount  of  impurity. 


A 


yOU-’r 


UvJll 


i;sRp  rhanap  iinder  nressure  of  E 


Figure  4  is  an  example  of  this  for  the  benzene- 
cyclohexane  system.  It  is  pressure  versus  fraction  of 
solid,  but  it  is  about  the  same  as  pressure  and  volume 
decrease  between  and  A2  .  The  curve  A  is  99.5% 
benzene.  The  curve  B  is  100%  benzene.  Clearly,  curve  A 
shows  the  transition  pressure  elevation  phenomenon. 


Fraction  Solid  (JC) 


Figure  4.  Transition  pressure  of  pure  (B)  and  impure 
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(A)  benzene  vs.  fraction  solid 


Common  Processes 


There  are  two  methods  to  apply  the  driving  force  to 
a  crystallization  process.  The  first  and  perhaps  most 
obvious  is  to  remove  heat  to  drop  the  temperature.  This 
is  the  more  obvious  due  to  the  frequency  of  occurrence. 

In  most  areas  ice  is  a  common  winter  sight.  The  second 
driving  force  is  to  increase  the  pressure.  It  should  be 
known  that  the  freezing  point  of  most  liquids  increases 
with  the  addition  of  pressure.  It  should  also  be 
realized  that  the  liquid  will  freeze  with  the  addition  of 
enough  pressure. 

Common  Temperature  Processes 

With  an  understanding  of  phase  diagrams,  some  of  the 
common  crystallization  processes  may  be  discussed.  The 
first  and  the  simplest  is  normal  freezing.  This  process 
slowly  freezes  the  impure  mixture  to  a  well  defined 
interface.  Freezing  proceeds  until  the  desired  amount  of 
product  is  frozen  and  the  remaining  mother  liquor  is 
discarded.  This  process  is  also  known  as  progressive 
freezing.  To  identify  a  progressive  freezing  process,  an 
interface  perpendicular  to  the  longest  dimension  of  the 
crystallizer  will  be  found. 


*  U"s«  -v.-^^"  »>-  v:  HJ  T  *  f  ■  L  W  w  >  if  ■.« 


1  ^  ’T. 
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Gouw^^  attempted  to  optimize  batch  normal  freezing 
by  using  programmed  temperature  with  time.  His  study 
found  that  impurity  deposit  was  not  affected  early  in  the 
freezing  cycle.  However,  late  in  the  cycle,  the 
impurities  in  the  deposited  solid  could  be  minimized  with 
the  programmed  temperature  approach.  Normal  freezing  has 
been  used  to  purify  continuous  solid  solutions. 

Fractional  melting  is  another  method.  It  uses  slow 
freezing  followed  by  slow  melting.  In  this  process  the 
impurity  freezes  last  and  melts  first.  The  impurity 
remains  in  concentrated  pockets  in  the  melt.  The  pockets 

drain  to  leave  a  progressively  purer  solid  behind. 

12 

Glassgow  and  Ross  cite  some  considerations  for  use  of  a 
fractional  melting  process.  First,  the  system  should 
segregate  simply  by  cooling.  No  solid  solutions  are 
formed.  Second,  the  system  mass  should  be  slowly  cooled 
and  reheated  close  to  equilibrium  conditions,  avoiding 
any  inclusion  of  the  impurity.  They  recommend  use  of  a 
closed  system  to  avoid  problems  involved  with  handling  of 
toxic,  odorous  or  unstable  compounds.  They  also 
recommend  an  adequate  driving-force  to  remove  the  mother 
liquor  adhering  to  the  crystals.  Desirable  systems  would 
produce  large  crystals  and  have  a  mother  liquor  of  low 
viscosity.  Aston  and  Mastrangelo^ ,  also  recognizing 
these  factors,  used  compressible  vanes  to  squeeze  the 
mother  liquor  from  the  crystals  and  in  doing  so,  reduced 
the  equilibrium  time. 
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The  first  industrial  crystallizer,  the  Proabd 

refiner,  was  originally  used  to  purify  napthene. 

Economics  of  an  industrial  process  are  very  important  and 

equilibrium  is  traded  for  productivity.  The  separation 

factor  obtained,  K,  will  be  less  than  that  which  could  be 

achieved  by  equilibrium.  For  supervision  and  labor 

requirements  to  be  minimized,  the  process  is  operated  in 

long  cycles  which  vary  from  24  -  48  hours.  Half  of  this 

time  is  used  for  cooling,  the  remaining  half  for  warming, 

draining  and  refilling.  The  crystallizer  is  a  large 

container  with  a  heat  exchange  network  capable  of  uniform 

heat  transfer.  The  flow  rate  of  the  heat  exchange  medium 

is  kept  high  to  maintain  a  small  temperature  difference, 

,5°  C,  and  uniform  conditions.^® 

The  product  liquid,  when  drained,  will  be  of  varying 

composition.  Many  fractions  could  be  removed,  however, 

for  practical  considerations  only  three  are  removed.  The 

mother  liquor  is  the  first  removed  and  is  later  discarded 

or  recycled  to  a  unit  process  before  the  crystallizer. 

The  second  is  a  crystallizer  recycle,  usually  of  purity 

equal  to,  or  greater  than  that  of  the  feed.  The  final 
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fraction  is  the  pure  product  fraction. 

Common  Pressure  Processes 

It  is  obvious  that  heat  transfer  problems  must  be 
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solved  to  operate  the  Proabd  refiner,  however,  this  paper 

is  concerned  with  the  pressure  driven,  adiabatic, 

crystallization  process.  There  are  two  methods  to  be 

considered  to  carry  out  the  crystallization.  The  first 

is  a  gas  pressurized  process.  The  idea  with  this  method 

is  to  use  the  gas  to  apply  the  work  and  crystallize  the 

liquid.  This  is  desirable  due  to  the  low  cost  of  the 

crystallizer.  There  would  be  no  working  parts  or  piston. 

The  use  of  a  reciprocating  compressor  would  introduce 

more  cost  but  perhaps  less  than  the  large  piston  of  the 

second  process.  There  is  no  known  process  which  is 

operated  in  this  manner.  The  second  method  is  hydraulic. 

A  large  piston  at  low  pressure  drives  a  smaller  piston  in 

the  crystallizer.  Working  in  this  manner,  M.  Moritoki 

and  co-workers  have  developed  a  pilot  plant  capable  of 
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producing  industrial  quantities.  Their  process  is 
adiabatically  pressurized  which  is  accompanied  by  a 
temperature  rise.  This  Kobe  Steel  group  found  that  the 
initial  state  and  increase  in  pressure  would  determine 
the  temperature,  solid  fraction  and  mother  liquor 

1 8 

concentration.  Figure  5  is  a  diagram  of  the  process. 

This  process  operates  in  five  steps.  Figure  6  is  a 
pressure-temperature  diagram  of  operation  of  the  process. 
First,  A,  the  piston  is  raised  to  introduce  the  feed 
followed  by  the  second  step,  B  -  C,  the  adiabatic 
pressurization.  Third,  C  -  D,  the  liquor  is  removed  by 
filtering  through  a  mesh.  The  fourth  step,  D  -  E,  is  the 


Figure  6.  Pressure-Temperature  relationship  for  Kobe 
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Steel  adiabatic  crystallization  process 


a-  equilibrium  line  of  pure  substance 
b-  equilibrium  line  for  initial  concentration 
c-  equilibrium  line  for  eutectic  concentration 


20 


r»_W  V'*. 
% 


k>  iT*.  w^l."*  :.^ 


Pressure  is  rather  uniform  throughout  the  crystallizer 
and  the  operating  cycle.  This  is  the  major  factor  which 
reduces  the  minimum  operating  cycle  time  and  uniform 
impurity  distribution. 
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Figure  9.  Vapor-liquid  phase  diagram  for  illustration 


.-v  .-^  ■■■^  P.' *y  »T^.vr».i  \i r.'-  "J  r.»-J  VJy.r:".  y.  y; wj»-.-r.’ r:i-;'»Ti:yr 


23 

As  illustrated  in  Figure  9,  high  increases  in  purity 
can  be  realized  by  distillation  in  the  middle  range  but 
as  an  increasingly  pure  liquid  is  fed  in  each  stage  a 
lower  increase  in  purity  is  realized-  An  increasing 
number  of  stages  are  required  to  obtain  a  product  more 
than  99%  puce.  The  result  is  a  much  higher  energy  cost. 


Minimum  Work 


In  an  economic  comparison,  energy  costs  are  also 
important.  The  minimum  energy  of  separation  related  to 
the  each  process  is  different.  The  minimum  work  of  each 
process  is  needed  to  compare.  First  distillation  is 
studied.  Then  the  minimum  work  applying  to  the  high 
pressure  process  is  examined. 


Minimum  Work  of  Distillation 

King^^  classifies  distillation  as  a  potentially 
reversible  process.  That  can  be  viewed  similarly  to  a 
heat  engine.  The  column  operates  by  taking  in  heat  and 
producing  a  separation.  Then,  the  reversible  heat 
entering  would  be: 


TqJ/To 


(7) 
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The  heat  leaving  would  be: 


(T, 


-  To)/Tl 


(8) 


Heat  entering  the  high  temperature  engine 
Heat  entering  the  low  temperature  engine 
High  temperature 
Low  temperature 
Reservoir  temperature 


The  net  work  would  be  the  difference  between  the  heat 
entering  and  the  heat  leaving  the  system. 


"n  ■  Oh  -  °L  I^L  - 


Assuming  no  mechanical  work,  the  energy  differences 
between  the  product  and  feed  will  be  small.  The  heat 
flows  reduce  to: 


«  Ik  ■  «  ^ 

'  «  •  •  •  *4  •  ■  ' 

*  -  *  -  ***  * 
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Qh  -  Ql  -  Q  (10) 

and  equation  (9)  simplifies  to; 

Wn  -  Q  Tq  (1/Tj^  -  1/Tjj)  (11) 

If  the  condenser  is  considered  the  low  temperature 
region  and  the  reboiler  considered  the  high  temperature 
region,  then  equation  (11)  becomes: 


"n  -  Oc  -^0  '1  ’'c  -  ‘ 


112) 


is  determined  by  jsinq  m.rimum  reflux 
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1  -  X, 


II  n 


(13) 


1 


Also,  wher. 


'  '  e  mi niraum  reflux 


A  is  the  latent  heat  of  vaporization. 


'‘■min 


Also  for  a  close  boiling  mixture  >  d. 

combining  (15)  and  (16)  leads  to: 


Q^,  -  /(«r  -  1)  F 


This  and 


(17) 


o<  and  Tv  can  be  determined  using  the  Clausius  - 

O 

Clapeyron  equation: 


d  In  P 


d  1/T 


(18) 


The  overhead  temperature,  T^,  for  a  component  in  a 
complete  separation  is  close  to  the  boiling  point  of  that 
component  at  the  column  pressure.  Likewise,  T^  is  close 
to  the  boiling  point  of  the  less  volatile  component.  The 
vapor  pressure  of  the  more  volatile  component  at  the 
bottoms  temperature  would  be  oC  times  the  column  pressure. 
By  integrating  the  Clapeyron  equation,  the  result  is: 


In 


?i/R  (1/T^  -1/Tj^) 


With  the  assumption  of  an  Or  close  to  1  ,  or  a  close 
boiling  mixture.  In  0<  equals  -  1.  The  o<  for  the 

benzene-cyclohexane  system  is  1.059  at  0°C,  1.027  at 
25°C,  and  1.003  at  50°C.  The  is  some  temperature 
dependence  but  the  assumption  of  close  to  1  remains 
valid . 


• 


•  A 
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The  assumption  of  and  by  combining  (11),  (17)  and 


(19)  : 


W„  -  R  F 
n  o 


For  any  reflux  above  minimum,  the  equation  would  be 
increased  by  a  factor  of  actual  vapor  rate  divided  by 
minimum  vapor  rate: 


W  -  (Actual  vapor  rate )/(Minimum  vapor  rate)  R  F  T 
-  ''■Actual  *  "’/"■min  *  4)  ■>  r 


Ninimum  Work  of  Pressurization 


The  minimum  work  for  the  pressurization  process  is 
much  simpler.  If  the  assumption  is  made  that  the 
compression  is  done  adiabatically ,  reversibly,  the  solid 
is  pure,  and  all  solid  is  recovered  as  pure  product,  then 
the  minimum  work  is: 


P  dv 


Given  an  equation  of  state  for  the  solution,  the  work  can 
easily  be  calculated. 


.V.v.v.v 
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IV.  Design  of  Process 


To  do  a  specific  comparison  of  distillation  to 

high  pressure  crystallization,  a  material  or  system 

must  be  chosen.  Benzene  of  90%  and  cyclohexane  of  10% 

were  used.  Benzene  is  an  important  building  block  for 

styrene  or  phenol  in  many  modern  synthetic  fiber  and 

19 

resin  processes.  Now  that  a  system  is  chosen,  the 
production  rate  must  be  determined.  Presently,  several 
conventional  benzene  production  plants  are  under 
construction  worldwide.  These  plants  are  designed  for 

C  O  1  g 

annual  production  rates  between  10  and  10  Kg/year. 

Also,  Chaty^  compares  various  industrial  processes 

where  the  annual  flow  rates  for  benzene  production  are 

10^-10^  Kg/year.®  With  these  data,  an  annual 

6 

production  rate  of  1x10  Kg/year  was  chosen. 

The  third  decision  is  to  determine  a  reasonable 

cycle  time  for  the  high  pressure  crystallization. 
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Mortiki  and  co-workers  indicate  a  minimum  cycle  time 
for  their  process  of  2  minutes,  allowing  10  seconds  for 
pressurization  and  equilibrium,  40  seconds  for 
filtration  and  the  remaining  time  for  sweating  and 
compacting . 


A'--.  -  v; -- 
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Correspondence  with  Hr.  Mortiki  suggests  a  time  of 
several  minutes.  The  absence  of  data,  and  the  decision 


to  use  an  inert  gas  for  the  compression,  required  a 
high  estimate  of  30  minutes  per  cycle. 

A  phase  diagram  for  this  system  is  also  necessary. 
Figure  10  is  the  phase  diagram  for  the  benzene- 
cyclohexane  system. 

Finally,  an  initial  temperature  of  25°  C  was 
selected  for  the  adiabatic  process. 


Crystallizer  Design  with  Gas  Pressurization 

A  basis  has  been  completed  so  the  process 

equipment  may  now  be  sized.  The  separation  of  the 

cyclohexane  is  first  determined.  Moritoki  published  a 

relationship  of  solid  phase  concentration  to  initial 
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concentration  for  the  impurity.  Figure  7. 

At  25°  C,  the  C„/C^  is  found  to  be  .027,  which 
s  o 

corresponds  to  a  product  of  99.7%  benzene. 

An  outlet  concentration  must  be  selected  for  the 
mother  liquor.  There  is  an  optimum  between  increase  in 
pressure  and  reduction  in  waste.  The  outlet 
concentration  was  selected  to  be  75%  benzene, 

2 

corresponding  to  a  maximum  pressure  of  1400  Kg/cm  ,  as 

1 8 

seen  in  Figure  10.  Assuming  8,000  hours  of  operation 


and  16,000  cycles  per  year,  the  product  rate  per  cycle 
is  62.5  Kg  or  137.75  lbs. 


Thus  far,  no  specifics  about  the  process  have  been 
discussed  and  the  process  layout  need  not  be  known. 

The  crystallizer  is  central  to  the  design.  The 
means  by  which  to  increase  the  pressure  must  also  be 
provided.  Therefore,  a  compressor,  followed  by  a  heat 
exchanger,  will  be  necessary.  The  crystallizer  feed 
stream  may  also  need  to  be  precooled  to  remove  the 
heats  of  compression  and  fusion.  Hence,  the  flow 


diagram  would  be  Figure  11.  The  crystallizer  is  shown 
in  Figure  12. 


REFRICERATION 


Figure  11.  Flow  diagram  for  gas  pressurized 


Mass  Balance 


Knowing  concentrations  in  all  streams  and  the 
product  flow  rate,  a  mass  balance  may  be  completed  on 
the  crystallizer  loop. 


Product  *  62.5  Kg/cycle 

Total  Mass:  Feed  «  Product  -f  Waste  (22) 

Cyclohexane  Mass:  C^- Feed  -  C^*  Product  +  Waste  (23) 

(.l)(Feed)  -  (.002)(62.5)  +  (.25)(Waste) 

by  solving  (22)  and  (23). 

Feed  -  102.9  Kg/cycle 
Waste  -  40.4  Kg/cycle 

Crystallizer  Volume  at  P^ 

At  atmospheric  pressure,  the  feed  volume  is  known 
and  the  system  volume  can  be  calculated. 

Volume  - 

Vg2  ■  1.138  cm^/gm 
V^jj  ■  1.284cm^/gm 

Volume  -  (102.9)(  .9)(1.138)  +  ( 102 . 9 ) ( . 1 ) ( 1 . 284  ) 
Volume  «  118.62  L 

Crystallizer  Volume  at  P 

g 

For  benzene,  Gehrig  and  Lentz  report  a  modified 


Tait  equation  for  P  in  MPa: 
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V  -  V  (B  +  P) 

-2 -  -  A  In  -  (25) 

V  (B  +  P  ) 


In  the  25°  C  range: 

A  -  .1845 
B  -  194.5  NPa. 


For  cyclohexane,  another  version  of  the  Tait 
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equation  is  reported  for  P  in  Kp/cm  : 


V  -  V  B  +  P 

2 -  .  C  log  -  (26) 

''o  ® 

In  the  25°  C  range: 

C  -  .2212  , 

B  -  72  3.3  Kp/cm“^  . 


P  is  1400  Kg/cm^  at  25°  C  and  75%  benzene.  Also, 

P^  is  1  atm  and  the  V^'s  are: 
o  o 


V 

V 


o  BZ 
o  CH 


1.138  cm?/gra 
1.284  cm'^/gm 


Each  equation  of  state  uses  different  units  of 


pressure. 


P  -1400  kg/cm^  -137.0  MPa  -1350  atm  -1395  kp/cm^ 


V  —  V 

'^O  BZ  1  BZ 


1  BZ 


1  BZ 


194.5  +  137.0 

.1845  In  - 

194.5  +  .101 

1.036  cm^/gm 


From  Bridgman^,  the  change  in  volume  of  solid 
benzene  A V/V^  is  .0902.^ 


A  V 


V 


o 


V  * 

S  BZ 

V  -  V, 

_0_CH _ _1_CH 

CH 

'^l  CH  - 


-  V 

-  .0902 

^o 

.943  cm^/gm 

920  +  1395 

-  .2272  In  - 

920 

1.167  cm^/gm 


At  the  maximum  pressure,  the  volume  will  be 
composed  of  the  product  amount  of  solid  and  the  liqu 
quantity  of  waste. 


Volume  -  solid  benzene  +  liquid  benzene 
+  liquid  cyclohexane 

Volume  -  (167.5) ( .943)  +  ( 40 . 4 ) ( 1 . 036 ) 

+  (40.4) ( .25) (1.167) 

-  102.1  L 

Liquid  Volume  -  43.2  L 

Nitrogen  Needed 

The  volume  decrease  between  P  and  can  now  be 

o 

calculated  and  the  amount  of  liquid  is  known.  The 
amount  of  nitrogen  needed  to  be  compressed  is 
calculable.  This  solubility  data  on  Nitrogen  in 
benzene  is  available  at  120  atm  and  25° 


It  is 


understood  that  the  operating  pressure  of  the  system  is 
much  higher.  But,  at  the  highest  pressure  the  system 
would  be  mostly  solid.  The  solubility  of  nitrogen  in 
the  solid  would  be  small.  These  data  will  be  assumed 
close  for  this  liquid  system. 

Xj,-  ■  .0551  mole  fraction 
-  .02050  mass  fraction 

So,  the  nitrogen  dissolved  is: 

mass  solution  «  (  .  02050 ) ( 40 . 4 ) 

-  .828  Kg 

AV  -  118.6  -  102.2 

-  16.5  L 

PAV  -  ZnRT 

-  40.3 
T^  -  2.36 

2  -  2.2 

(1350) (16.51)  -  (2.2)  n  (82.06)(298) 

n  ■  .4140  Kg  mole 
-  11.59  Kg 

Total  mass  nitrogen  =  12.68  Kg 

Due  to  the  assumed  solubility  this  may  be  a  high 
estimate . 


The  volume  (SCF)  of  gas  can  also  be  calculated; 

( 1)  V  -  (  .4530) (  .7306)  (637) 

-  463  SCF 


Pressure  Vessel  Design 


The  crystallizer,  a  high  pressure  vessel,  has  a 
large  number  of  regulations  involved  with  its  design. 
The  American  Society  of  Mechanical  Engineers,  (ASME), 
have  established  rules  and  standards  governing  safety, 
design,  fabrication  and  inspection.^  The  desired 
vessel  is  a  cylindrical  shell  with  a  50  cm  ID  and  2:1 
ellipsoidal  heads.  The  maximum  pressure  is  19,700  psi 
and  material  of  construction  is  A-302B  with  an 
allowable  stress  of  26,500  psi.^  All  joints  will  be 
fully  radiographed  in  accordance  with  UW-ll(a), 
Reference  1.  The  joint  efficiency  is  1.0  when  joints 
are  fully  radiographed. 

Shell  circumferential  stress  is  calculated  for 
thick  cylindrical  shells  from  Appendix  l,l-2(a)(l), 
Reference  1. 


■5  -  U 

(29) 

E  +  P 

(30) 

E  -  P 

'(26500  +  ISTOO)] 

1 

,(26500  -  19700)J 

40.2  cm 


Shell  longitudinal  stress  is  also  calculated  for 


thick  cylindrical  shells  in  accordance  with  Appendix  1, 
l-2(a)(2).  Reference  1. 

t  -  R  (Z-^  -1)  (31) 

Z  -  (P/(S  X  E)  +1)  (32) 

t  -  (25) ( (19700/26500  +1)*^  -1) 

«  8.01  cm 

The  largest  of  the  two  values  is  used.  The  shell 
thickness  would  be  40.2  cm.  The  heads  are  to  be  2:1 
ellipsoidal  in  accordance  with  UG-32(d),  Reference  1. 

P  X  D 

t  - -  (33) 

2  S  X  E  -  .2  P 

(19700) (50) 

(2) (26500) (l)-( .2)(19700) 

-  20.07  cm 

To  proceed  further  in  sizing  the  crystallizer  the 
length  of  the  vessel  must  be  determined.  This  is  found 
by  subtracting  the  volume  if  the  heads  from  the  total 
volume  needed  and  dividing  that  difference  by  the 
crystallizer  cross-sectional  area.  Once  the  length  is 
known,  an  approximate  volume  of  the  construction 
material  of  can  be  found.  Since  pressure  vessels  are 
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roughly  priced  by  total  weight,  these  calculations  will 
complete  the  cost. 

Vg  -  4/3  <r  a  b  c  (34) 

The  variables  a,  b,  and  c  are  the  lengths  of  the  semi- 
axies . 

V_  -  (4/3)'rr(25)  (25)  (12.5)/(1000) 

32.7  1 

V-  -  135.3  -  (2)(32.7) 

-  69.9  1 

-  (69900)/  (25)^ 

-  35.6  cm 

The  volume  of  building  material  of  the  shell  and 
the  heads  must  be  determined.  The  volumes  calculated 
using  the  outside  radius  are  subtracted  from  the  volume 
inside  radius. 

Vpjs  *  V  L  ( (R+tg)  -  R)  (35) 

-  35.6  ((25+  40.2)  -  25) 

-  4.053  X  10^  cm^ 

V„H  -  (4/3)ir  (  (a+tj^)  (b+tjj)  (c+tjj)-abc)  (36) 

-  (4/3)  fr  ((45.1)(45.1)(32.6)-25^  12.5) 

-  2.778  X  10^  cm^ 

Vtotal  ■  ’-SO®  * 

-  .29  Ib/in^,  , 

-  1.77  X  10"^  Ib/cm'’ 

Mass  -  1.204  X  10^  lb 


Applying  data  from  Figure  6,  Reference  6,  the  cost 
per  pound  would  be  $  1.20  and  an  additional  $  .60  for 
stainless  steel  construction. 


Total  cost  ■  $  21,700 


Minimum  Work 


Finally,  the  minimum  work  for  the  compression  of 
the  liquid  mixture  must  be  calculated.  This 
calculation,  as  stated  earlier,  is  simply  the  integral 
of  PdV.  The  work  is  calculated  in  3  parts. 

The  integrals  of  PdV  for  liquid  benzene  and  liquid 
cyclohexane  are  summed  over  each  range  of  specific 
volume.  The  integral  of  PdV  is  also  calculated  for  the 
transition  of  benzene  from  liquid  to  solid  and  the 
compression  of  the  solid.  The  following  equations  are 
used; 


for  liquid  benzene! 


V.  -  V 


A  In 


B  +  P 


B  +  P. 


for  liquid  cyclohexane: 


V  -  V 


B  +  P 


B 


The  benzene  liquid  -  solid  transition  is  known.  Data 
from  Bridgman  for  solid  benzene  was  developed  into  the 
line  for  the  P-V  relationship  over  the  interval  680  <  P 
<  1500  kp/cm^.^ 

P  -  101,500  (Vq  -  V)/Vq  -  15,520  (37) 

For  illustration,  large  step  sizes  will  be  used. 
More  reasonable  steps  are  used  with  a  design  program  to 
be  used  with  Flowtran,  only  a  few  points  are  reported 
and  graphed.  Data  was  only  available  for  one  point 
between  the  pressure  of  initial  crystallization  and  the 
highest  pressure  achieved.  These  data  give  the  line 
above . 
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N  Table  1 

K  Results  of  Work  Done  on  Liquid  Components 


^  Pressure 

( Kq/cm^ ) 

Pressure 

Benzene 

Cyclohexane 

average 

AV 

Pave  A  V 

AV 

Pave  AV 

. 

( atm) 

( L/Kg ) 
10"^ 

(atm  L/Kg) 

(L/«g 

10"^ 

)  (atm  L/Kg) 

■ 

'a 

1  1 

b 

h 

P  100 

49.0  1 

.000 

.5 

1.273 

.6 

1 

200 

145.6 

.920 

1.3 

1.152 

1.7 

300 

242.7 

.855 

2.1 

1.068 

2.6 

400 

339.8 

.805 

2.7 

.967 

3.3 

500 

436.9 

.745 

3.3 

.908 

4.0 

1  600 

534.0 

.700 

3.7 

.845 

4.5 

1  700 

631.1 

.655 

4.2 

.790 

5.0 

i 

Total  work  per  kilogram 

17.8 

r_  (atm 

L/Kg) 

t  800 

728.2  1 

.88 

13.7 

.748 

5.4 

1  900 

825.2 

.66 

5.4 

.707 

5.8 

1 

'  1000 

922.3 

.63 

5.8 

.673 

6.2 

1100 

1019.4 

.61 

6.2 

.629 

6.4 

1200 

1116.5 

.58 

6.5 

.605 

6.7 

;  1300 

1213.6 

.55 

6.7 

.573 

7.0 

1350 

1286.4 

.27 

3.5 

.280 

3.6 

P  Total  work  per  pound 

65.6 

62.9 

1 

(atm 

V 

L/Kg) 

Table 

2 

S 

a  Results  of  Work  Done 

on  Solid  Benzene 

i 

r 

Pressure 

k 

^  2 

Pressure 
( average ) 

TRANS 

> 

<1 

> 

Pave  AV 

PaveAVqiQ.  Mc 

S  (Kg/cm  ) 

% 

( atm) 

( L/Kg ) 

( atm 

L7Kg7“ 

W  680 

1596 

1.045 

■  1400 

1003 

1794 

1.051  .00494  4 

.95 

179.9 

F 

F  2500 

• 

1992 

1.058 

& 

Total  work  per  pound  184.9 

(atm  L/Kg) 


These  totals  can  be  multiplied  by  the  quantity  of 


each  component  of  the  mixture  to  give  an  estimate  of 
the  work  required  for  the  compression. 


Work  «  Solid  Benzene  +  Liquid  Benzene 

+  Liquid  Cyclohexane  (38) 

-  (184.9)(62.5)  +  (65.6)(30.2)  +  (17.8)(62.5) 

+  (62.9)(10.3) 

«  15,300  L  atm 

-  1500  BTU 

-  10.25  BTU/Kg 


Tables  1  and  2  give  reasonable  estimates  for  the 
work  required.  A  more  accurate  model  was  constructed 
by  writing  a  FORTRAN  program  to  do  the  above 
calculations.  Table  3  gives  the  result  of  that  program 
and  Figure  13  graphically  shows  the  pressure- 
temperature  relationship. 
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Table  3 

Work  on  Benzene-Cyclohexane  System 

Pressure(atm)  Work(BTU)  Temperature 

Cycle  lb  Product  Dif ference(°F) 


23.6 

7 

.1 

.1 

158.9 

75 

.6 

.8 

294.2 

174 

1.3 

1.9 

429.5 

286 

2.1 

3.2 

564.9 

405 

2.9 

4.5 

700.2 

524 

3.8 

5.8 

835.5 

645 

4.7 

7.1 

970.9 

765 

5.6 

8.5 

1016.0 

862 

6.3 

9.5 

1061.1 

961 

7.0 

10.6 

1106.2 

1044 

7.6 

11.6 

1151.3 

1116 

8.1 

12.4 

1196.4 

1180 

8.6 

13.1 

1241.5 

1238 

9.0 

13.7 

1286.6 

1291 

9.4 

14.3 

1331.7 

1407 

10.2 

16.7 

Precooler 


From  the  data  it  is  clear  that  the  heat  of 
compression  (work)  and  heat  of  fusion  must  be  removed. 
This  will  allow  the  system  to  crystallize  at  the 
desired  temperature  of  25  °C.  A  cooler  must  be 
designed  to  remove  1407  BTU  per  operating  cycle. 
Again,  an  overall  heat  transfer  coefficient  of  26 
BTU/lbmol  ft^  °F  will  be  used.  The  temperature 
difference  will  be  20  °F.  This  will  be  carried  out  in 
a  stirred  jacketed  tank. 

This  process  needs  some  means  of  cooling  the  feed 


stream  by  removing  1407  BTU  in  2/3  of  a  cycle.  A 
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HIGH  PRESSURE  CRYSTALLIZER 

BLOCK  -  BHPC  NAHE  -  CRYS 

OPERATING  PRESSURE  >  13S4.292ATH 

CRYSTALLIZER  UOLUHE  *  117.779L 

SPECIFIC  VOLUME  OF  BENZENE  Ccc/ill 

LIQUID  AT  PO  1.13800  LIQUID  AT  P  1.05327 

SOLID  AT  P  0.95638 

SPECIFIC  VOLUME  OF  N-HEXANE  Ccc/il 

LIQUID  AT  PO  »  1.28400  LIQUID  AT  P  >  1.16857 

SHELL  THICKNESSCIN3»1S.S5492SHELL  MASSCLB3«  8693.15 

HEAD  THICKNESS! IN3-  7.83211HEAD  HASSCLBT-  3954.87 

CRYSTALLIZER  COST  -  $  22766.44 

MINIMUM  WORK  PER  HOUR  *  2815.598BTU 

MINIMUM  WORK  PER  POUND  PRODUCT  9*483 

NITROGEN  DISTRIBUTION  CLB3 

IN  SOLUTION  s  6.60IN  GAS  PHASE  -  23.37 

CONCENTRATION  OF  IMPURITY  (N-HEXANE) 

IN  »  0.097414PR0DUCT  *=  0. 0016S6M0THER  LIQUOR  =  0.750000 


Figure  14.  Flowtran  output  for  gas  pressurized 
system 
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refrigeration  plant  of  .40  ton-day  capacity  would  be 
sufficient.  Figure  13-81,  Reference  30,  has  data  on 
NH^  refrigeration  plants.  The  ammonia  expansion  is 
assumed  at  constant  temperature  so  the  assumed 
temperature  difference  is  tolerable.  A  system  capable 
of  10  ton-day  would  cost  $  17,000. 

Refrigeration  Plant  , 

Cost  -  17,000  ( .40/10) 

-  $  2,500 

Heat  Transfer 

A  -  A  H/U  AT  (  39) 

-  (1407)(6)/(26)(20) 

-  15.8  ft*^ 

Cost  for  heat  transfer  area^®  -  $1700 

The  volume  of  the  tan)c  should  equal  that  of  the 
crystallizer,  40  gallons. 

Cost  mixed  tank^^  «  $4000  for  100  gallons 

Cost  -  4000  ( 40/100) 

-  $2300 

Total  Precooler  Unit  Cost  »  $  6,500 


The  first  law  simplifies  to: 


-  Wn  -  (H2  -  P2‘'''2^  “  ^”1  "*‘l‘'^l^  "  Am  (43) 

-  -  (6316-(150) (51.88) )  -  ( ( -38 )-( 1 . 5 ) ( 1659 ) ) 

"  -  (38)(453) 

-  -  4  48.3  BTU/cycle 

The  minimum  work  the  compressor  must  do  is  the  sum 
of  the  minimum  work  on  the  liquid  and  the  minimum  work 
on  the  nitrogen  gas.  This  sum  multiplied  by  a  72% 
effiency  for  the  compressor  and  divided  by  the  time  the 
compressor  operated  would  give  an  approximate  power 
requirement  for  the  compressor  motor. 


-  W  -  Work  on  Gas  +  work  on  Liquid  (44) 

"  -  448.3  +  1407 

-  1855  BTU 

Power  -  1855/( .72) (600) 

-  6.07  hp 


2 

Also  with  similar  data  on  nitrogen,  the  outlet 
temperature  can  be  calculated. 


AH 


c; 

AtP 

AT 


(H,  -  H. )/E 
C  '‘dT 

6316  J/gmol 
38  J/gmol 
35.0  J/gmol-  K 


(6316-38)/( 32.0) ( .72) 
195  _C 


382  °F 


(45) 

(46) 


The  compressor  is  sized  by  a  correlation  of  cost 
vs  power. The  required  power  is  lower  than  any 
graphed,  so  a  cost  size  extrapolation  factor  is  used. 


rJTT 


For  100  hp,  the  cost  of  the  unit  is  $30,000  at 


5000  psi 


This  must  be  extrapolated  to  8.1  hp  and 


20,000  psi  discharge. 

Cost  -  30,000  (6.1/100)°-®  (20,000/5,000)°-° 

-  $  9,700 

The  heat  exchanger  may  be  designed  knowing  the 
compressor  design,  the  heat  capacity  of  nitrogen  and 
the  outlet  temperature  of  the  gas.  The  compressed  gas 
temperature  must  be  lowered  to  the  crystallizer 
temperature,  298°  K. 


Mass  flow  rate  -  .4140  Kgmol/10  min 
-  .04140  Kgmol/min 

,H  -  H,  -  H, 

-  898  -  6316 
»  -5418  J/s 
•  .09930  BTU/s 


The  heat  flow  rate  must  equal  the  change  in  energy 
in  the  nitrogen  stream. 


ah  -  U  A  A  T, 


U  for  gases  and  water  is  in  the  range  of  2-50 
BTU/lbmol • f t^- °F , ^°  an  average  of  26  BTU/lbmol • f t^- °F 
will  be  used  to  calculate  an  approximate  area.  The 
coolant  will  be  water  inlet  65°  F  and  outlet  90°F. 


AT,  -  (  AT.  -  AT^)/ln  (At./AT^) 

-  ( 12.0i292)/in  (12.0/29i)  ° 

-  87.7  °F 


fk  .-mm 


.  "A" 

V.-.'.'V*.-*.'.-.' 


---A*  -V- 

.  <  .r,  .dR.V-  . 


s.  ••  V 

.  -.-•'a'Cn 

I**-*.'.*. 


•.*  ■ 


'  a'*  «>  **• 
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A  «  AH/(UAT,  )  (39) 

-  ( .0993Q)(36oO)(6)/(26)(87.7) 

-  .94  ft^ 

With  area  and  pressure  known,  a  corresponding  cost 
can  be  determined  form  graphical  data.  Such  data  is 
available  in  Figure  14-17,  Reference  30.  With  an  area 
of  1.0  ft^  and  stainless  steel  construction,  the  cost 
is  $1300.00.  An  adjustment  factor  for  higher  pressure 
is  needed.  Using  the  size  adjustment  factor  of  .44  for 
pressure  vessels  and  since  the  change  in  pressure  is  a 
change  in  size,^^  this  will  be  assumed  to  apply: 

Cost  -  1300  (  20,000/600) 

-  $  6000 

Nitrogen  Storage 

A  storage  vessel  must  be  used  to  recover  the 
nitrogen  let  down  from  the  crystallizer.  The  gas  will 
be  stored  at  200  psi  and  enough  gas  for  5  cycles  will 
be  stored.  Some  loss  will  result  from  gas  in  solution 
and  the  volume  of  gas  left  at  200  psi. 

Total  -  .414  Kgmol  x  5 
*  2.07  Kgmol 

(13.3)(V)  -  (2.07)  (87.06) (298) 

-  3,805  L 

.  134  ft^ 


V 
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This  tank  will  be  designed  like  the  crystallizer. 
A  radius  of  36  inches  and  height  of  60  inches  will  be 
used  for  a  finished  volume  of  141  ft^. 


Shell  circumferential  stress: 

t  -  P  X  R/(S  X  E  -  .6P) 

-  (200)(36)/((15000)(l)  -  (.6)(200)) 

-  .484  in 

Shell  longitudinal  stress: 

t  »  P  X  R/{S  X  E  -  .4P) 

-  (200)(36)/( (2)(1500)(1)  -  (.4)(200)) 

-  .239  in 

2  heads  elliptical  heads  10% 

t  -  P  X  D/(2S  X  E  -  .2P) 

-  (200)(72)/( (2)(15000)(1)  -  (.2)(200)) 
■  .480  in 


The  mass  of  the  tank  is  necessary  to  complete  the 
cost  determination. 


Volume  shell  =  4y'((R+t)^  -,R^)  L.,  (35) 

=  «r  (  (  36.484)^  -  36^)  (60) 

=  6610  in^ 

Volume  head  »  (4/3)  if*  ( a+t )  ( b+t )  (  c+t )  -  abc  (36) 
»  (4/3)  4r  (36.24)(18.24)-  (36)(18) 

=  5300  in'* 

Density  *  .284  Ib/in^ 

3 

Total  volume  ■  17,200  in 
Total  mass  »  4,890  lb 


Cost  from  Figure  6,  Reference  7,  is  $0.80  per 
pound  for  low  alloy  steel. 


(40) 


(41) 


(33) 


Cost  -  $  3,900 


•y-yy 

'-*>-' 
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Monsanto's  Flowtran  has  an  extensive  library  of 

32 

unit  operations  to  model  most  processes.  Often  there 
is  a  need  to  develop  new  programs  for  unit  processes 
not  in  the  library.  A  program  to  be  used  in 
conjunction  with  Flowtran  has  been  developed.  This 
program  does  the  work  of  the  previous  design.  This  is 
a  batch  process  though  and  as  such  does  not  convert 
well  with  Flowtran.  This  creates  problems  with  using 
the  program.  A  listing  of  the  Flowtran  program  can  be 
found  in  the  Appendix.  Figure  14  is  the  output  of  the 
new  Flowtran  block.  As  seen  the  program  reports  many 
of  the  intermediate  values  as  well  as  the  cost  of  the 
vessel , 

Distillation  Column 

The  distillation  column  for  comparison  with  the 
crystallizer  can  also  be  modeled  by  Flowtran.  The 
overall  rates  will  be  the  same  but  they  must  be 
converted  to  continuous  rates.  The  parameters  of  the 
column  were  first  determined  by  a  preliminary  Flowtran 
model  which  used  the  Underwood  equation.  The  parameters 
were  reflux  ratio,  number  of  trays,  feed  plate,  bottoms 
pressure,  condenser  pressure,  and  distillate  to  bottoms 
ratio. Those  parameters  were  used  with  a  more  rigorous 
model.  The  rigorous  model  selected  the  diameter.  The 
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TITLE  DISTILLATION  COLUHN  FOR  COHPARISON 
PROPS  21222 
r<ETR  CH  DZ 
PRINT  POUNDS 

BLOCK  DIST  DISTL  SI  S3  S2 
PARAM  DIST  1  26*9  60  4S  15  13  *3982  000 
BLOCK  DICO  CDSTL  DIST 
PARAH  DICO  1000 
LBS  SI  1  45*3  408*2 
TEHP  SI  70 
PRESS  SI  14 
END  CASE 
.END  JOB 


DISTILLATION  COLUHN  FOR  COHPARISUN 

PHYSICAL  PROPERTY  OPTIONS 
ANTOINE  VAPOR  PRESSURE 
REDLICH-KWONC  VAPOR  FUGACITY 
CORRECTED  LIQUID  FUGACITY 
SCATCHARD-HILDEBRAND  ACTIVITY  COEF 


Figure  16.  Input  to  Flowtran  simulation 
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DISTILLATION  COLUHN  FOfI  CQHFARISQN 


OIST  -  OISTL  -  TOTAL  CONOCNSCC: 

FEED  •  SI  •  iOTTONS  •  S3  >  OVERHEAD  •  S3 

reflux  RATIO*  28. to.  NO.  OF  PLATES  •  AO..  FEEO  PLATE*  AS. 

FRAC  OVHD  •  0.3t8.  FEED  FRAC  VAPOR*  0.000 

condenser  DUTY*  O.ttAOEFOADTU/MR.  TEMP*  176.S3F>  PRES*  IS.OOPSIA 
REDOILER  DUTY*  0.tl33CFOADTU/HR.  TEMP*  169. 73F.  PRES*  IS.OOPSIA 


DICO  (DISTO  COST  FOR  UNIT  DIST 

nuhder  or  diameters  i. 
NUHDER  of  trays  60.0 
HOC  FACTOR  FOR  TRAYS  1.00 
LOADING  AT  TOP  O.BS 
PLATE  SPACING  AT  TOP  18. 
SURFACE  TENSION  AT  TOP  30. 
FOAMING  FACTOR  AT  TOP  1.00 
OVERALL  U  FOR  CONDENSER  100.0 
COOLANT  TEMP  RISC  20.00 
COOLANT  HEAT  CAPACITY  1.00 
CONDENSATE  SUBCOOLINO  0.00 
TYPE  OF  CONOENSER  -1. 
HOC  factor  for  CONDENSER  1.00 
COOLANT  COST-C/MOL  3.00 
HEATING  FLUID  MAX  TEMP  700.00 
MIN  TEMP  OIF  FORFLUX  1.00 
TYPE  OF  REDOILER  -1. 
PRESSURE  FACTOR-RCDOILER  1.00 
REDOILER  UPDATE  FACTOR  124.0 
STRESS  IN  COLUMN  SHELL  14000. 
EXTRA  UCIGHT  ITEMS  2300. 


REFLUX  RATIO  20.900 
TYPL  OF  tray  3. 
TRAY  COST  UPDATE  FACTOR  124.0 
LOADING  AT  tOTTOH  O.OS 
PLATE  SPACING  AT  BOTTOM  18. 
SURFACE  TENSION  AT  BOTTOM  20. 
FOAMING  FACTOR  AT  BOTTOM  1.00 
COOLANT  TEMP  IN  80.00 
MINIMUM  TEMP  APPROACH  10.00 
coolant  SPECIFIC  GRAVITY  1.00 
PRESSURE  factor-condenser  1.00 
TUBE  LENGTH  FACTOR-COND  1.00 
CONDENSER  UPDATE  FACTOR.  124.0 
HEATING  FLUID  COST-C/MLB  56.00 


REDOILER  HEAT  FLUX  12000. 
HEATING  FLUID  DELTA  H  1000. 
HOC  FACTOR  FOR  REBOILER  1.00 
TUBE  LENGTH  FACTOR-REDL  1.00 
SKIRT  HEIGHT  10.0 
CORROSION  ALLOUANCE  0.000 
shell  COSTrtFLD  O.SO 


CONDENSER - 

CONDENSER  DUTY.  DTU/HR  *  0.894037SSEF06 

LOG  MEAN  TEMP  DIFFERENCE  >  86.13  DEO  F 

AREA  *  103.80  80  FT 

COOLANT  FLOW  RATE  *  89.404  GPM 

CONDENSER  COST  *  «  1260.00 

coolant  COST  *  6/MR  0.16 

REBOILER  - 

reboiler  DUTY.  BTU/HR  -  0.9132692SEF06 

AREA  *  76.11  SO  FT 

steam  FLON  rate  *  913.269  LD/HR 

REBOILER  COST  •  6  1030.00 

STEAN  COST  *  6/HR  0.51 


COLUMN  DIAMETERS-  COLUHN  HEIGHT.  FT  *  100.00 


BELOU  FEED  2.0000  SHELL  THICKNESS.  IN  *  1.6250 

ABOVE  FEED  3.0000  SHELL  THICKNESS.  IN  *  1.6250 

COST  OF  COLUMN  SHELL  *  6  22270.00 

TRAY  COST  *  6  3620.00 


TOTAL  UTILITY  COST*  6/HR  0.67 

COLUMN  CAPITAL  COST  *  6  28200.00 


Figure  17. 


Flowtran  output  for  distillation  column 


56 


57 


desired  result  was  a  column  optimized  between  operating 
cost  and  capital  cost  while  producing  a  product  of 
equal  purity  to  the  crystallization  process.  The 
process  diagram  is  shown  in  Figure  15.  The  Flowtran 
input  is  in  Figure  16.  The  results  of  the  Flowtran 
design  are  shown  in  Figure  17. 

Crystallizer  Design  with  Hydraulic  Compression 

An  alternate  to  the  high  pressure  design  can  also 
be  proposed.  Figure  18  is  a  diagram  of  the  alternate 
design.  It  is  still  a  high  pressure  process  but, 
instead  of  an  inert  gas  pressurizing  the  system  the 
pressure  is  supplied  hydraulically  and  intensified. 

The  crystallizer,  Figure  19,  is  a  double  ended  piston 
and  cylinder  arrangement.  On  one  end  an  incompressible 
fluid  is  pumped  to  supply  the  force  on  the  piston.  The 
force  is  equal  to  pressure-area  of  piston  product.  The 
area  of  the  crystallizer  piston  is  much  smaller  than 
the  fluid  piston  allowing  a  large  increase  in  pressure. 

P2  -  Pi  (Aj/A2)  (42) 

The  design  of  these  two  pressure  vessels  and  the  piston 
will  be  similar  to  the  pressure  vessels  presented 
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Crystallizer  Similarities 

This  process  has  several  similarities  to  the 
previous  crystallization  process.  First/  the  mass 
balances  and  cycle  times  will  be  the  same.  This  will 
also  force  the  crystallizer  volume  to  be  the  same.  The 
minimum  work  on  the  liquid  will  also  be  equal  due  to 
the  above  reasons.  The  minimum  work  will  necessarily 
be  less  than  the  gas  compression  process,  since  there 
is  no  work  on  gas  compression.  This  will  also  make  the 
operation  expense  of  this  process  lower  than  the  first 
process.  The  final  similarity  is  the  precooler.  The 
work  done  on  the  fluid  and  the  heat  of  fusion  liberated 
will  be  the  same  so  the  precooler  size  and  cost  will  be 
the  same . 


Pressure  Vessel  Design 

The  material  of  construction  will  also  be  the  same 
as  the  first  design,  302-SS.  The  low  pressure  fluid 
will  be  200  psi.  The  area  difference  of  the  two  piston 
heads  can  be  calculated  with  the  vessel  diameters  and 
wall  thicknesses  of  the  shells  and  heads.  Subscripts  1 
and  2  are  for  the  high  and  low  pressure  ends, 
respectively . 
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Area 

20,000  -  200  (A,/A,) 

A-  -  320  cm*'  i 
A“  -  32,000  cvT 

Radius 

r-  ■  10.1  cm 
ri  *  101  cm 
»  370  cm 

Thickness  and  volume  of  high  pressure  vessel  are 
calculated  as  the  previous  high  pressure  vessel,  p.37- 
40. 


16.2  cm 

XT  (26.3  -  10.1)  370 
685,000  cm-* 


t„,  -  8.10  cm  - 

VIJ2  -  (4/3)  ft  (  (18.2)  ^(13.2) 

-  (10.1)2(7.6)1 

-  4500  cm^ 


Using  the  area  and  the  volume,  the  crystallirer 
dimensions  can  be  calculated.  The  volume  of  solid  is 
the  total  volume  with  the  liquid  volume  subtracted.  The 
volume  of  solid  will  allow  calculation  of  the  piston 
length,  since  it  is  the  maximum  travel  of  the  piston. 
The  volume  of  the  low  pressure  vessel  may  also  be 
calculated  from  this  length. 


solid 


58.9  L 


-  (58.9)(1000)/(320) 
piston  _  200  cm 

Thickness  and  volume  of  low  pressure  shell 
V, 


SHI 


(200X  32,000) 
-  6400  L 


'SHI 
t. 


'HI 


1.36  cm 
.673  cm 


SHI 


fr  (102.36^  -  101^)  200 


173,500  cm' 


HI 


(4/3)  ((101.673)^(51.173) 


-  (101)‘‘(50.5)  ) 
.3 


■  58,00  cm' 

Piston  Volume  -  (200)  ir  (10.1)^  +  (10)  4t  (101)^ 
Total  Volume  Material 


1.306  X  10®  cm^ 


Material  density  ■  1.77  x  10~^  Ib/cm^ 


Material  Cost  -  $  1.80 
Total  Cost  *  $  41,600 


Hydraulic  Pump 


The  pump  for  the  low  pressure  fluid  should  be 
specified  to  pump  the  total  amount  of  fluid  in  2/3  of 
the  cycle  time.  It  can  be  sized  from  graphical  data. 
Figure  13-42,  Reference  20. 


Flow  rate  -  ( 6400 )/( 20 ) ( 3 . 178 ) 
-  22  gpm  @  200  psi 


psi 


The  cost  found  is  $  2000  for  a  discharge  of  100 
This  needs *to  be  adjusted  to  200  psi. 

.6 


Cost 


2000  (200/100) 

-  $  3100 

The  pump  also  needs  a  10  hp  motor  which  costs  $ 


;  ■  %  •  «» 
%  *  .  -  V  • 

/-'A 


-■•V  V 


•.i'wvV.'.i'.-vJ 


650  from  Figure  13-56,  Reference  30. 
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Three  processes  have  been  designed  to  perform  the 
same  separation.  Two  designs  crystallize  the  material 
as  means  of  separation.  The  third  is  a  distillation 
process,  which  separates  by  vaporization.  The  summary 
of  major  equipment  and  operating  costs  for  these 
processes  is  located  in  Table  4.  The  summary  of 
operating  cost  is  in  Table  5. 


I 


! 


Table  4 

Summary  of  Crystallization  and  Distillation  Processes 


1.  Crystall 

izer  with  gas  compression 

Unit 

Specifications 

Cost( 1985 

Crystallizer 

302-stainless  steel 

Shell  r  »  25  cm 

1  ■  35.6  cm 
t  »  40.2  cm 
Ellipsoidal  heads 
r.  «  25  cm 
■  25  cm 
r|  ■  12.5  cm 
t  ■  20.1  cm 

53,200 

Compressor 

6.1  hp,  20,000  psi 
discharge, 200  psi  feed 

29,500 

Heat  Exchanger 

7.0  ft^.  Water  coolant 
in  at  65°  F,  T  «  25°  F 

8,500 

Gas  Storage 

Carbon  Steel 

Shell  r  »  36  in 

1  *  60  in 
t  -  .484  in 
Ellipsoidal  Heads 
r  =  36  in 
r  *  36  in 
r  =  18  in 
t  =  .240  in 

6,600 

Precooler 

NH,  Plant-  .40  ton-day 
Mixed  Tank,  40  gal 

Heat  Exchanger 

3,200 

3,300 

2,400 

Pi 


Table  4,  Continued 


2.  Crystallizer 

with  Hydraulic  Compression 

unit 

Specifications 

Cost( 1985 

Crystallizer 

302-stainless  steel 
Hydraulyc  Intensifier 
and  crystallizer  vessel 
200  psi  r  »  101  cm 

1  -  200  cm 
t  -  .693  cm 

20,000  psi 

r  =  10.1  cm 

1  =  cm 

. ’^HEAD  1 

^cuiTTf  16.2  cm 

Piston®?^s!^101  cm 

1  -  10  cm 
r  -  10.1  cm 

1  »  200  cm 

Hydraulic  fluid, 1600  gal 

105,500 

8,000 

Pump 

200  psi  discharge 

22  gpm 

4,300 

Motor 

10  hp 

930 

Precooler 

NH,  Plant-  .52  ton-day 
Mixed  Tank,  40  gal 

Heat  Exchanger 

Pump 

3,200 

3,300 

2,400 

815 

Total  Equipment 

Cost 

$"12  4;  9  4 '5 

3.  Distillation 

Column  for  Comparison 

unit 

Specifications 

Cost( 1985 

Condenser 

Reboiler 

Trays 

Column 

103.8  ft^ 

76.1  ft^ 

60  sieve  trays 
with  downcomers 

100  ft  high 

18  in  between  trays 

3,500 

2,800 

9,800 

60,300 

Total  Equipment 

Cost 

$  76,400 

-w’  J_V sJ. JV 


Table  5 

Summary  o£  Anual  Operating  Costs 


1.  Crystallizer 

with  gas  compression 

Cost 

160 

16 

606 

4,540 

Utility 

Refrigeration 
Cooling  Water 
Electric 
( Compressor ) 
Nitrogen  gas 

Ammount 

.52  ton-day/day 

30  gal/cycle 

6.1  hp  for  10  min/cycle 

46.0  SCF/cycle 

Total  Utilities 

Cost 

$  5,325 

2.  Crystallizer 

with  Hydraulic  Compression 

Utility 

Ammount 

Cost 

Refrigeration 

.52  ton-day/day 

160 

Electric  (Pump) 

10  hp  for  10  min/cycle 

999 

Total  Utilities 

Cost 

$  1,159 

3.  Distillation 

Column 

Utility 

Ammount 

Cost 

Steam 

914  Ib/hr 

10,500 

Cooling  Water 

89.4  gpm 

2,200 

Total  Utilities  Cost  $  12,70t) 


To  compare  the  different  processes,  there  must  be 
several  assumptions.  First,  it  must  be  assumed  that  this 
process  will  be  a  part  of  a  plant,  which  produces  benzene  and 
further  processes  a  pure  benzene  product.  An  example  could 
be  ring  closure  of  n-hexane  producing  cyclohexane. 
Dehydrogenation  of  the  cyclohexane  produces  the  benzene. 
Benzene  can  be  further  processed  to  a  monomer  such  as  phenol 
or  styrene. 

The  crude  benzene  is  first  purified  by  a  crude 
distillation  to  90%  purity.  The  benzene  is  then  further 
purified  by  one  of  the  designed  processes.  For  this  reason, 
land,  yards,  service  facilities,  and  working  capital  have 
already  been  provided.  The  focus  is  economic  evaluation  of 
the  purification  necessary  to  the  plant  operation. 
Installation  must  be  added  regardless  of  the  design  selected. 
It  will  be  calculated  from  a  30%  factor  of  the  purchased 
equipment  cost.  The  percentage  factors  of  purchased 
equipment  cost  are  referred  to  as  Lange  Factors. There  are 
Lange  factors  to  determine  expenses  such  as  instrumentation, 
electrical,  and  piping.  It  is  assumed  that  these  expenses 
would  be  higher  for  distillation  than  for  high  pressure 
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crystallization.  Also  the  Lange  factors  reported  are  for  all 
chemical  processing  industries.  The  factors  would  not 
compensate  for  the  less  expensive  instrumentation  of  the 
crystallization  and  for  this  reason  are  not  included.  The 
addition  of  these  services  would  increase  each  process 
proportionally  which  would  probably  not  be  the  best  cost  to 
compare . 

The  present  values  of  the  installed  equipment  and  the 
present  value  of  10  years  of  operating  costs  are  computed  and 
the  sum  compared  to  select  the  most  economical  process. 

Table  6 

Installed  Equipment  Cost  and  Operating  Expense 

Process  Equipment  Installation  Operation  Total 

Annual  10  years 

1  108,015  32,400  5,325  40,380  180,795 

2  124,945  37,500  1,159  8,790  171,300 

3  76,400  22,900  12,700  96,900  196,200 

4  40,900  12,300  2,400  18,200  73,400 

Process  1,  Crystallizer  with  gas  pressurization 

Process  2,  Crystallizer  with  hydraulic  pressurization 

Process  3,  Standard  distillation  column 

Process  4,  Proabd  refiner 

The  total  cost  for  the  hydraulically  pressurized 
process  is  the  lowest.  This  would  probably  be  the  process 
chosen.  There  may  be  other  factors  involved  beside  the  cost. 
The  plant  may  use  many  standard  pumps  and  may  not  use 
compressors.  This  might  influence  the  decision  to  go  with 


the  more  expensive  compressor  design.  There  may  also  be 
reliablility  data  on  compressors  and  piston  systems  that  may 
make  the  use  of  one  more  attractive.  The  development  of  a 
pilot  plant  for  each  of  the  crystallizer  processes  could 
answer  many  of  the  questions  about  reliablility. 


Table  7 

Separation  Costs  of  Crystallizer  and  Distillation  Designs 

Process  Total  Feed  Separation  Cost 

Installed  Cost  Composition  ($/Kg) 

1  140,400  90%  .00533 

2  162,500  90%  .00116 

3  99,300  90%  .01270 

4  205,500 

(53,200)  98%  .00240 


Process  1,  Crystallizer  with  gas  pressurization 
Process  2,  Crystallizer  with  hydraulic  pressurization 
Process  3,  Standard  distillation  column 
Process  4,  Proabd  refiner 


Table  7  shows  the  designed  processes  and  a 
refrigeration  type  industrial  process.^  The  Proabd  refiner 
as  discussed  earlier  operates  on  batch  cycle  and  uses 
programmed  freezing  and  warming  to  separate  the  feed.  As 
shown  in  Table  7,  the  comparison  of  costs  for  separation 
clearly  show  that  the  high  pressure  processes  have  a  sizable 
economic  advantage  over  distillation  and  crystallization  by 


refrigeration.  The  refrigeration  process  was  reported  for 
economic  comparison  for  complete  plants.  The  installed  cost 
was  estimated  with  Lange  factors  and  is  reported  in 
parenthesis  in  the  Table  7.  The  feed  composition  is  also 
important.  The  Proabd  refiner  is  the  least  costly  but,  is 
fed  a  98%  pure  mixture.  The  effect  of  raising  the  feed 
impurity  concentration  would  greatly  increase  the  capital 
investment  and  raise  the  operating  costs.  The  crystallizer 
would  need  a  larger  design  to  freeze  the  increased  mother 
liquor  and  this  would  also  force  the  operating  cost  higher, 
probably  higher  than  the  high  pressure  processes. 

Before  implementation  of  the  high  pressure 
crystallization  process,  several  other  questions  need  to  be 
studied.  The  first  is  to  gather  data  on  specific  systems  of 
interest.  A  multipurpose  pilot  plant  should  be  developed  to 
study  various  systems  and  techniques  of  operation  and 
control.  If  the  gas  pressurization  design  is  selected,  there 
is  an  excellent  chance  of  a  patent  on  the  process.  The 
hydraulic  process  may  also  be  eligible  for  a  patent.  Kobe 
Steel  does  not  have  United  States  patents  on  the  process  they 
have  designed  so  it  may  be  possible  to  alter  or  improve  on 
their  design.  The  cycle  time  assumed  in  the  crystallizer 
design  was  30  minutes.  The  Kobe  Steel  research  concluded 
that  cycle  times  could  be  as  low  as  2  minutes.  They  report 


. ' . . . » .  • 


.-.'l-’.-.'vV/ 


|l'.  /|p- 

‘v*  •..* 


si 


A/*  • 

*/  ,v  •», 


1^ 


■  -'I  '**.  - 


*  -*•  a 

.  ^  *  a  a  .  “ 


.>  .*•  . 


71 


the  equilibrium  to  be  established  in  as  low  as  10  seconds. 

If  the  cycle  time  were  reduced,  the  present  design  could  be 
reduced  in  size  and  much  less  costly  to  construct.  The 
construction  of  a  pilot  plant  could  confirm  this  as  well  as 
determine  where  various  separation  processes  could  be 
improved  upon. 


VI .Conclusions 

The  results  of  this  study  of  high  pressure 
crystallization  can  be  summarized  as  follows: 

1.  The  hydraulic  compression  crystallization  process 
is  the  least  expensive  of  the  processes  study.  It  is  capable 
of  separating  the  feed  benzene  stream  for  a  capital  cost  of 

$  162,445  and  operate  for  $  1,159/year.  This  is  the  lowest  of 
the  processes.  Table  4. 

2.  The  hydraulic  pressurization  may  have  some  patent 
problems.  The  Kobe  Steel  research  group  has  a  process  which 
is  similar.  They  have  an  operating  pilot  plant  and  several 
patents  granted  and  pending. 

3.  The  gas  pressurization  process  is  close  to 
competitiveness  with  the  hydraulic  process.  It  is  lower  in 
capital  cost  but,  is  higher  in  operating  costs.  The 
operating  costs  could  be  reduced  by  using  a  diaphragm  to 
allow  compression  by  filtered  air  instead  of  nitrogen. 
Reduction  of  operating  costs  could  also  be  realized  from  a 
more  complete  recovery  of  the  nitrogen. 

4.  The  use  of  distillation  for  this  separation  is  not 
recommended.  The  high  cost  required  to  operate  the  process 
would  be  a  great  deal  more  expensive  for  the  expected 
equipment  service  life. 


VII. Recommendations 

The  future  of  this  project  should  proceed  as  follows; 

1.  Prior  to  implementation  of  the  high  pressure 
crystallization  process  data  on  specific  systems  should  be 
gathered.  A  multipurpose  pilot  plant  should  be  built  to  study 
these  systems  and  to  develop  techniques  of  operation  and 
control.  If  the  gas  pressurization  design  is  selected,  there 
is  an  excellent  chance  of  a  patent  on  the  process.  The 
hydraulic  process  may  also  be  unique  enough  for  a  patent. 

Kobe  Steel  does  not  have  United  States  patents  on  the  process 
they  have  designed  so  it  may  be  possible  to  alter  or  improve 
on  their  design.  The  cycle  time  assumed  in  the  crystallizer 
design  was  30  minutes.  The  Kobe  Steel  research  concluded 
that  cycle  times  could  be  as  low  as  2  minutes.  They  report 
the  equilibrium  to  be  established  in  as  low  as  10  seconds. 

If  the  cycle  time  were  reduced  the  present  design  could  be 
much  less  costly  to  construct.  The  construction  of  a  pilot 
plant  could  confirm  this  as  well  as  determine  where  various 
separation  processes  could  be  improved  upon. 

2.  Other  improvements  could  include  an  optimum  switch 
from  distillation  to  crystallization.  The  crystallizer  will 
produce  the  high  purity  at  low  cost  but,  only  when  feed  a 
high  purity  material.  There  could  be  an  optimum 
concentration  to  switch  separation  techniques  from 


distillation  to  crystallization.  In  the  case  of  this  study, 
the  optimum  is  probably  at  a  feed  of  higher  than  90%  pure. 
With  a  more  pure  feed  the  crystallizer  could  be  designed  at  a 
lower  pressure  and  smaller  cycle  capacity.  The  pressure 
reduction  would  result  from  good  productivity  at  a  lower 
ultimate  mother  liquor  concentration.  The  crystallizer 
capacity  may  also  be  reduced  because  with  a  more  pure  feed 
the  equilibrium  in  the  system  will  be  achieved  more  quickly. 
The  result  of  these  factors  would  allow  a  less  expensive 
design  for  the  crystallizer  as  well  as  lower  operating  costs 
due  to  the  reduction  of  mother  liquor. 

This  work  has  possibly  led  to  a  patentable  process. 
The  application  for  patent  has  begun.  If  granted  it  would  be 
well  advised  to  fund  a  pilot  plant  which  could  answer  many 
questions  about  process  parameters  (i.e.  cycle  time,  feed 
concentrations,  product  purity)  having  influence  on  the 
overall  cost.  The  pilot  plant  could  be  used  to  gather  data 
on  many  systems  and  to  determine  optimal  concentrations  to 
switch  from  distillation  to  crystallization.  The  pilot  plant 
could  also  answer  some  question  about  the  reliability  of  the 
equipment  used  in  the  crystallization. 


VIII .BIBLIOGRAPHY 


1. 


2. 


3. 


4. 


5. 


6. 


7. 


8. 


9. 


American  Society  of  Mechanical  Engineers. 

Section  VIII,  Rules  for  Construction 
of  Pressure  Vessels,  Division  1. 

New  York:  United  Engineering  Center, 
1983. 

Angus,  S.,  de  Reuck,  K.H.,  and  Armstrong,  B. 

International  Thermodynamic  Tables  of 
the  Fluid  State-  6  Nitrogen.  New  York; 
Pergamon  Press,  1977 . 

Aston,  J.G.,  and  Mastrangelo,  S.V.R.  "Improved 
fractional  apparatus."  Analytical 
Chemistry,  April  1954,  26(4),  p.  764- 


Bowden,  W.W.  Classical  Thermodynamics  for 

Engineers.  Terre  Haute,  Indiana :  Rose- 
Hulman  Institute  of  Technology,  1981. 

Bridgman,  P.W.  Proceedings  of  the  American 
Academy  of  Arts  and  Sciences,  1949 , 
it;  p.‘*129.  - 

Chaty,  J.C.,  and  Wilcox,  W.R.  "Economics  of 

fractional  solidification."  In  M.  Zief 
and  W.R.  Wilcox  (Eds.)  Fractional 
Solidification  (Vol.  1).  New  York: 
Marcel  Dekker,  Inc.,  1967. 

Clark,  F.D.,  and  Terni,  S.P.  "Thick  walled 
pressure  vessels."  In  Modern  Cost 
Engineering;  Methods  and  Data.  New 
York :  Me  Graw-Hill  Publishing  Co . , 

1979,  p.  339-342. 

Denbigh,  K.  The  Principles  of  Chemical 

Equilibrium  (4th  ed.).  New  York ; 
Cambridge  University  Press,  1983, 
p.  197-200. 

Gelrig,  M.,  and  Lentz,  H.  "PVT  data  for  benzene 
in  the  range  5  to  300  MPa  and  323  to 
683  K . "  Journal  of  Chemical 
Thermodynamics,  9.  1977 .  57  4  4  5-4  50. 


Gel'perin,  N.I.,  Laphenkov,  G.I., 

Nosov,  G.A.,  and  Taran,  A.L. 
"investigation  of  the  process  of 
crystallization  of  a  melt  on  surfaces 
with  the  aid  of  resistance  networks." 
International  Chemical  Engineering, 
l8(ll),  January  1978,  p.  Il5-ll8 . 

Gerard,  W.  Solubility  of  Gases  and  Liquids;  A 
Graphical  Approach.  New  York;  Plenum 
Press ,  1976^  p.  2^4 . 


Glassgow,  A.R.,  Jr.,  and  Ross,  G.  "Purification 
of  substances  by  a  process  of  freezing 
and  fractional  melting  under 
equilibrium  conditions."  Journal  of 
Research  of  the  National  Bureau  of 
Standards ,  September  1958,  ( Research 
Paper  1113) ,  ^(3),  p.  137-142. 

Gouw,  T.H.  "Programed-temperature  normal 

freezing."  Separation  Science,  1967, 
2(4),  p.  431-417. 

Guthrie,  K.M.  "Estimating  the  cost  of  high 

pressure  equipment."  In  H.  Popper  (Ed.) 
Modern  Cost  Estimating  Techniques. 

New  York:  Me  Graw-Hill  Publishing  Co. , 
1970,  p.  147-150. 

Hoffman,  H.C.  "Pressure  vessel  primer." 

Machine  Design,  52(23),  9  October  1980, 
p.  168-171.  — 

"Hydrocarbon  processing  industry  boxscore." 

Hydrocarbon  Processing,  Section  2, 
February  1984,  p"!  14-43 . 

King,  C.J.  Separation  Processes.  New  York; 

Me  Graw-Hill  Publishing  Co.,  1971. 

Kobe  Steel,  LTD.,  New  Materials  Division. 
"Technical  presentation  on  high 
pressure  crystallization  process." 

Leyerrle,  J.D.  "Supply/demand  of  benzene  in 

1990."  Chemical  Engineering  Progress, 
80(5),  5  May  1984,  p.  l3-l6, 

Molinari,  J.G.D.  "The  Proabd  refiner."  In  M. 

Zief  and  W.R.  Wilcox  (Eds.)  Fractional 
Solidification  (Vol.  1).  New  York: 
Marcel  Dekker,  Inc.,  1967. 


Molinari,  J.G.D.  "The  Newton-Chambers  process." 
In  M.  Zief  and  W.R.  Wilcox  (Eds.) 
Fractional  Solidification  (Vol.  1). 

New  York:  Marcel  Dekker,  Inc.,  1967. 

Moritoki,  M.  "A  method  of  fractional 

crystallization  under  high  pressure." 
International  Chemical  Engineering, 
^(3),  July  l9S0,  p.  394-401. 

Moritoki,  M.,  and  Fujikawa,  T.  "Liquid-solid 
phase  diagram  of  p-  and  m-cresol 
mixtures  under  pressure."  In  S.J.  Janie 
and  E.J.  de  Jong  (Eds.),  Industrial 
Crystallization  84.  Amsterdam, 
Netherlands :  Elsevier  Science,  1984, 
p.  369-372. 

Moritoki,  M.  "Crystallization  and  sweating  of 
p-cresol  by  application  of  high 
pressure."  In  S.j.  Janic  and  E.J. 
de  Jong  (Eds.),  Industrial 
Crystallization  84.  Amsterdam, 
Netherlands :  Elsevier  Science,  1984, 
p.  373-376. 

Moritoki,  M. ,  Kitagagawa,  K.,  Onoe,  K.,  and 

Kaneko,  K.  "  A  pilot  study  of  p-cresol 
crystallization  by  adiabatic 
application  of  high  pressure."  In  S.J. 
Janic  and  E.J.  de  Jong  (Eds.), 
Industrial  Crystallization  84. 
Amsterdam,  Netherlands :  Elsevier 
Science,  1984,  p.  377-380. 

Moritoki,  M."  Phase  diagram  sent  in  personal 
correspondence."  19  August  1985. 

Mulet,  A.,  and  Corripio,  A.B.  "Estimate  costs 
of  pressure  vessels  via  correlations." 
Chemical  Engineering  (New  York), 

88(20),  5  October  1981,  p.  145- 
rJ6, 148,150. 

Mullin,  J.W.  Crystallization (2nd  ed.). 

Cleveland ,  Ohio :  CRC  Press  Inc.,  1972. 

Osugi ,  J.,  Shimizu,  K.,  Yasunami ,  K.,  Moritoki, 
M.,  and  Onodera,  A.  "Liquid-solid 
transition  at  high  pressure  III." 

Review  of  Physical  Chemistry  of 
Japan"  38(2),  1968,  p.  90 . 


Peters,  M.S.,and  Timmerhaus,  K.D.  Plant  Design 
and  Economics  for  Chemical  Engineers 
(3rd  ed.).  New  York:  Me  Graw-Hill 
Publishing  Co.,  1960. 

Richman,  D.,  Wynne,  E.A.,  and  Rosi,  F.D. 

"Progressive  freezing."  In  M.  Zief  and 
W.R.  Wilcox  (Eds.)  Fractional 
Solidification  ( Vol .  1 ) .  New  York : 
Marcel  Dekker,  Inc.,  1967. 

Seader,  J.O.,  Seider,  W.D.,  and  Pauls,  A.C. 

Flowtran  Simulation-  An  Introduction 
( 2ncl  ed.).  Cambridge,  Massachusetts: 
CACHE,  1977. 

Taslimi,  M.  "The  PVT  data  of  20  organic  liquids 
and  their  relationship  to  molecular 
structure."  Doctoral  Thesis, 

Department  of  Mechanical  Engineering  at 
Technical  University  of  Hannover,  1969. 

Weast,  R.C.,  (Ed.).  Handbook  of  Chemistry  and 
Physics  (59th  ed. ) .  Boca  Raton, 

Florida :  CRC  Press,  Inc.,  1979. 

Wilcox,  W.R.,  and  Zief,  M.(Eds.).  Fractional 
Solidification  (Vol.  1).  New  York: 
Marcel  Dekker,  Inc.,  1967. 

Wilcox,  W.R.  "Introduction."  In  M.  Zief  and 
W.R.  Wilcox  (Eds.)  Fractional 
Solidification  (Vol"^  1 ) .  New  York : 
Marcel  Dekker, Inc.  ,  1967. 


Wolten,  G.M.,  and  Wilcox,  W.R.  "Phase 

diagrams."  In  M.  Zief  and  W.R.  Wilcox 
(Eds.)  Fractional  Solidification 
(Vol.  IT^  New  York:  Marcel  Dekker, 


APPENDIX 


c  This  is  a  fortran  subroutine  which  is 

c  designed  to  run  with  Flowtran.  tt  is 

c  block  BOdeling  a  batch  high  pressure 

c  crystallizer.  The  block  list  is  g  characters 

c  long, 

c  1  Block 

c  2  Naae 

c  3  BHPC 

c  4  Feed 

:  c  5  Gas  in 

(  c  6  product 

;  c  7  nother  liquor 

c  8  Gas  out 

c  The  paraaeter  list  has  7  pieces  of  data 

c  to  input. 

'  e  1  Taaperature 

c  2  Cycle  tiao 

c  On  Material  of  Construction 

c  3  naziaun  allowable  pressure(ata) 

c  4  Density  (g/cc) 

c  S  Joint  eficiency 

c  6  Cost/unit  weight  (5/lb) 

c  7  Radius  of  vessel  (in) 

c  The  retention  vector  contains: 

c  1  Operating  pressure 

c  2  Crystallizer  voluae 

c  specific  voluaes  of: 

c  3  Benzene  liqud 

c  4  Benzene  solid 

c  5  Cyclohexane  liquid 

c  6  Shell  thickness 

c  7  Head  thickness 

c  a  Shell  wieght 

c  9  Head  weight 

c  10  Crystallizer  cost 

c  11  Hiniaua  wock/hour 

e  12  Nitrogen  in  solutuion 

c  13  Nitrogen  in  space 

e  14  Concentration  in  feed 

c  15  Concentration  in  product 

c  The  units  throughout  the  prograa  ace  unifora. 

c  Thicknesses  or  lengths  are  in  inches, 

c  Weights  are  in  pounds, 

c  Density  is  aeasuted  in  g/cc 

c  Specific  voluaes  ace  in  cc/g. 

c  Pressure  is  in  ata, except  when  soae  equations 

c  rf  stats  ace  used.  The  units  are  stated  and 

c  converted  back  after  calculation, 

c  Tiae  IS  in  hours, 

c  All  flow  cates  are  in  pound  aoles. 

c 

SUBROUTINE  BHPC ( E 1 , G1 , Pi .W1 ,G2 , PAN , NET) 

DIHEMSION  rl<  31) ,Pl( 31),W1(31) ,G1( 31) ,G2( 31) ,PAN( 31) ,RET( 31) 
COHNON/B/  rN(3l),PR;31),ANLH(31),CI(3l),GO(31) 

COHNON/C/  P ,  VCPO  ,  VCP ,  VOBZ .  VLBZ ,  VSBZ ,  V0C6  ,  Vt,C6  ,  TSB , N5H ,  POI .  PRN 
COBNON/O/  TB , WB ,  WONRLB ,  CRCOST , WONKCT ,  WN2SOL ,  WN2SPC ,  CO ,  CS ,  Clll, 
COHHON/E/  T , C YCTH , S , BOU , E , WCOST , B 
DO  111  I«l,31 
FRd)-  ri(i) 

PHd)-  Pl(I) 

AMLR(l)-  Wl(I) 

GKD-  Gl(I) 


m-M.€6  718 


UNCLRSSIFIED 


ECONOniC  COnPARISOM  OF  A  HI8H  PRESSURE  CRVSTAUIZATION 
PROCESS  HITH  A  DISTILLATION  PROCESS(U>  ARHV  AILITARV 
PERSONNEL  CENTER  ALEXANDRIA  VA  J  E  DIETZ  87  JAN  8C 


M’CRnrnr 


G0(I)>02(II 

111  COHTIHUB 
NSTEr-  30 

ir(IIOUT.BQ.2)CO  TO  SO 
CrCTM-  OAK  3) 

T>  rA«(l) 

A-  OAJKS) 

S>  OAK  7) 

E-  PAJI(») 

CML  >  PAAUS) 

ROU'PAKll) 

»COST«  PAA(13I 
CALL  PRESS 
CALL  SEP 
CALL  COMPS 
CALL  CRWOl 
CALL  CRTS! 

CALL  MORR 
CALL  N2RB0 

SO  WRITE  I  6, 100) 

WRITEtO.lOl) 

WRITE) 6, 102)  PR(20) ,PR(20) .AMLR(20) 

WRITB(6,103)  GI(20), 00(20) 

WRITe(S,104)  T 
WRITE(6,10S)  CYCTM 
WRITEX.IOS)  S 
WRITE) 6, 107)  ROU 
WRITEX.IOS)  E 
WRITE) 6, 109)  WCOST 
WR1TE(0,110)  R 
WriTEX.IOO) 

WRITE(«,101) 

WRITE) 6, 120)  P 
WRirE(«,l21)  VCPO 
Wr;TE(6,13S) 

WR:Tf(6,122)  V08S,VLEZ 
WRITE) 0,123)  VSBS 
WRITE(0,130) 

WRITE (0,124)  VOCO.VLCO 
WR1TE)0,12S)  TSH.WSE 
WRITE)0,12«)  TB.WR 
WRITE) 0,127)  CRCOST 
WRITE(0,120)  WORRCT 
WRITCtO.lSQ)  MORELS 
WRITE)0,137) 

WRITC(0,129)  WN2S0L,WN2SPC 
WRITE(0,130) 

WRITB(0,130I  C0,CS,CML 

100  rORflAT('l'.lSE,'EI<!E  PRESSURE  CRTSTALLIZBR' ) 

101  PORNATdOS.’BLOCR  '  BBPC' , SX, ’NAHE  -  CRTS’) 

102  rORMATdOX, 'PESO  •  ’  ,A4, ’PRODUCT  •  ’,A4, 'HOTBER  LIQUOR  -  ’  ,A4 

103  FORMAT)  lOX, 'NITROOBM  OAS  IN  •  ’  ,A4, 'NITROOBM  GAS  OUT  •  ’,A4) 

104  rORMATdOX, 'OPERATING  TEMPERATURE  -  ',P0.2,'DBC  K' ) 

105  rORHATdOX, 'CYCLE  TIRE  >  '  ,PS. 4. 'B0UR8' > 

100  rORMATdOX, 'MAXIMUM  ALLOMIABLE  PRESSURE  -  '  ,  PIO.  4,  ’  ATM' ) 

107  PORMATdOX, 'MATERIAL  DENSITY  •  ’  .PE.S, '9/CC’ > 

100  rORMATi lOX, 'JOINT  BPPIBNCV  •  ',rS.S) 

109  PORMATdOS.'COST  PER  WBIGBT  •  '  ,P0.S. ’$/L»' ) 

110  rORMATdOX, 'RADIUS  OP  VESSEL  •  ’.PE.S,*  IN’) 

120  PORMATdOX, 'OPERATING  PRESSURE  •  ’  .PO.  3 ,  ’  ATM’ ) 

121  PORMATdOX,' CM YSTALLISER  VOLUME  -  '.PS.3,'L') 

122  PORMATdOX, 'LIQUID  AT  PO' ,  PO.  5,  lOX, 'LIQUID  AT  P’.PS.S) 


'  >^-t  K  *1.^ , 


f 


I 


•  '>■:  y.'?.'  .'■  '.•  J.  .'.'  i-.  y.  y,  f.- 


,'  IP.  'r.'wj"! 


123  F0IUIAT(16X.' SOLID  AT  P'.FS.S) 

12«  FOMATdSX,' LIQUID  AT  PO  >  PS .  S.lOX, 'LIQUID  AT  P  -  ',FS.S> 

125  POXHATdSX,' SHELL  THtCENESSdNl*' .PS.S, 'SMELL  MASS! LSI-' ,  P14 . 2 ) 

126  POMIATdex.'SBAD  THICKMBSSdHl-' .PS.S. 'HEAD  HASSt  LB)-' ,  P14 .2 ) 

127  rOMMATdSX, 'CMYSTALLtZER  COST  -  $',ri0.2) 

128  rORHATdSX. 'MINIMUM  NOBS  PEE  HOUR  -  '  ,P12 . 1, 'BTU' ) 

150  rORMATdSX. 'MINIMUM  MORE  PER  POUND  PRODUCT' ,  PS. 3) 

12P  PORMATdOX,  '  IN  SOLUTION  •  '.PS.  2. 'IN  GAS  PHASE  -  ',P8.2) 

130  PORHATdSX.'lN  -  '. PO . 6 PRODUCT  •  PS.S. 'MOTHER  LIQUOR  -  '.PS. SI 

135  PORMAT<  ISX. 'SPEClPtC  VOLUME  OP  BENZENE  |ce/9l'> 

)  13S  PORMATdSX.'SPECIPIC  VOLUME  OP  N-HEXANE  Icc/Ql' ) 

'  137  PORHATdSX. 'NITROdEM  DISTRIBUTION  (LBl'I 

130  PORMATdOX, 'CONCENTRATION  OP  IHPURITT  (N-BEXANE)') 

END 

C 

C 

C 

c 

SUBROUTINE  PRESS 

C 

c  This  subroutln*  calculate*  the  tiansltion 

c  pccasuca  for  banian*  given  a  taapacatur* 

c  and  uses  as  a  point  to  define  a  line 

c  knowing  the  slop*  of  the  phase  diagtaa. 

c  The  line  ia  Chen  used  to  deteraine  the 

c  highest  pressure  needed  to  provide  the 

c  given  outlet  concentration. 

C  FROM  P.94  REP. 55 

c  P-Po-  a( (T/To)*»c-ll 

c  and 

e  P  «  -3567  Cout  ♦  1 

c  Cout  -  outlet  concentration 

c  POI  •  intercept  dateralned  froa  pressure  (X  benxan*  -  1.0) 

e  [PI-  bars 

c  (Pol-  bar* 

c  a-  1.6  E3  bars 

c  c-  2. SO 

c  To-  noraal  bioling  point  at  Po 

c 

c 

COMMON/C/  P ,  VCPO . VCP . VOBZ . VLBI , VSBZ . VOCS , VLC6 , TSB ,N8R . POI . PEN 
CONHON/D/  TH , NN , NORRLB . CRCOST , NORKCY ,NM2SOL . NN2SPC, CO , C'  , CNL 
COMHON/E/  T.CYCTM.S.ROU.E.NCOST.R 
COHHON/B/  PR(  31 )  ,  PR(  31 ) .  ANLR(  31 )  ,CI  (  31 )  .GO(  31 ) 

A-3600 

C-2.S0 

T0-27fl.8 

PO-1. 01325 

TEMP-  T/TO 

PA-1»A«( (TEMP«-C)-1) 

PP-  PA/1. 01 325 
POI-  -2430  -  PP 
P  -  -2430-CML  -  POI 
RETURN 
END 
c 


c 

SUBROUTINE  SEP 
e 

c  Thia  subroutine  detaralnes  the  seperatlen  based 


fc' 

I 

•» 
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•*. 
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on  dots  fron  N.  Hotitoki,  *A  nothod  of  fractional 
crystlllation  undor  hlqn  piaasuca,*  Intetnatlonal 
Chaaical  Enqlnaarlnq.  Vol.20.  No. 3,  p.39T. 

Tha  Ca/Co  va.  praaauta  cutva  wat  appcoxiaatad 
by  a  atraight  lino  ovor  cho  900  to  000  kg/cal 
ranga: 

Ct/Co  -  .131«1  -  t.l74St-«  (»). 


CCMNON/C/  P . VCPO . vet ,V0BZ,VLiZ.VS»2, VOCS ,VIC6. TSH . WSM , PO I , PPM 
COlWON/ Q/  TH.im , NOPKia . CKCOS r .MORKCV .HM2SOL, IIM2SPC , CO , CS , CHI. 
COHHON/ e/  T , CY CTH , S , HOU , B , NCOST , ■ 

COHNON/Hy'  PR(  31 ) .  PR(11 1  ,>HLR(  31 )  .GK  31 )  .GO(  31 ) 

C0-rR( 3)*84/(PR(l)*7«*rR( 3)*S4) 

IP(P. 61.1000)00  TO  210 

CSC0«  .13181  -  ((1.17488-4)*Pi 

CSCO  •  .017 

CS*CSC0*C0 

RETURN 

END 


SUBRSUTtNE  COUPS 

This  aubcoutina  aitigna  tha  tlcw  rataa  of  tha  prpduct 
and  by  prodjct  atraaaa  by  uaing  aasa  balances. 

COMNON/C/  P,  VCPO. VCP.V0B:,Vl,PZ,VSBl,V0C6.Vt.C6, TSH, HiH, POI  .PRN 
COHHON/0/  TH,NH,NORKtB,CRCOST,NORRCY.NN2SOL,HN2SPC,CO,CS,CH: 
COMMON/ E/  T.CYCTM.S, ROU , C , NCOST , R 
CCMMON/B/  fR( 3l I .PRi 31  I ,AMLR( 31 » .Gl, 31 l,CO( 31 ' 

PROP  -  (IPR(1)*PR(3))*(1-CML  -PRC  3 ) )/) ( 1-CML) -CS i 

pr;i)  -  (i-:s'«pROP 

PRC  3)  •  CS'PROP 

PRINT*,  'PROP ‘.PROP 

AMLRCn  .  PRCD-PRCl) 

AM1,RC  3  )  •  PRI  3)-PR(  3  I 

RETURN 

END 


SUBROUTINE  CRYVOL 

This  aubroutina  calculacaa  tha  cryatilliaar 
voluaa  given  an  operating  ptasaura  and  tha 
physical  property  data  on  the  syatea  coaponenta. 

The  voluaa  ia  calculated  uaing  tha  aaauaption  that 
cyclohexane  and  bensene  are  independent  in  aaking  up  tha 
voiuRa. 


COM,MON/C/  P ,  VCPO ,  VCP ,  VOBl ,  VLBI ,  VSBE ,  VOCB  .VICE  ,  TSH  ,MSH ,  POI ,  PRN 
COMHON/D/  TH , NH , WORRLB . CRCOST , WORRCY ,NN2SOt , MNZSPC , CO , CS ,  CMC 
CONMON/E/  T, CYCTN, S, HOU, E, NCOST, R 
COMHON/B/  ERC  31 ) , PRC  31 ) , AMLRC  31 ) ,GI  (  31 )  ,GOC  31 » 

PO-l 

VOBZ*  1.138 


u  u  u  u 
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VOC«-  1.2S4 

AJISC«  •  CYCTN*ril(  l)*a4 

AMSBZ  •  cycTn*rt(  V)*'>a 

CALL  ESBZP 

CALL  escep 

VCPO-  (AMSBZ*V0aZ*AflSC4*V0ca>*.4S4 

vcp  -  (VLaz*AHLa( i)*7a*vsaz*pii(i)*T8*vLC6*ra( })*a4)*cycTii*.4S4 

RETUIW 

END 


suaaouTiHE  Esazp 
e 

c  This  •ubcoutin*  caXeulat**  th«  voIum  of  bonsona 

c  givan  a  praasura  and  taaparatura.  It  uaaa  publiabad  data  trom 

c  nanrad  Gahrlg  and  Harco  Lantt  in  Journal  of  Chaaleal 

c  Therucdynaaics,  1977,  9,  p.4S5-4S0,  and  P.W.  arldgaan 

c  tn  procaadingi  of  tha  Aaarican  Acadaay  of  Arts  and  Sclancaa, 

c  1949,  77,  p. 129-144.  Givan  a  praasura  abova  tha  transition 

c  praasura  tha  prograa  atsuaaa  ESt  solid  and  1S%  liquid 

c  bacausa  of  tina  lioitations  <  aquilibriua  is  not 

t  econouical  to  achiava  )  and  tapuritias  prasant  in  tha 

c  solution.  Tha  aquations  usad  ara: 

c 

c  {  Vo-V  )/V  •  Ain  (  (  B*P  )/(  B*P0  )  I 

c  IBI-  HPa 

c  tPl-  npa 

c  p.  449  Gahrig  and  Lantx 

c  At  fraaxtng  25  oC  and  480  kg.'ca2 

c  OaltaV/  V«  .0434  and  .1594 

e  p.  14C  Btidgaan 

e 

COnnON/C/  P  ,  VCPO  ,  vcp  ,  VOBZ  .  VLBZ  ,  VSBZ  ,  V0C4 ,  VLC4  ,  T5H ,  MSB ,  PO t ,  PEN 

ccn.9c/«/o/  IN , tm , worelb .  cacost , hobecy, i*n2soi , whzspc  ,  co  ,  cs . cbl 

COnnON/E/  T, CYCTH , S . BOO . E . MCOST , B 

connoN  /§/  rB(  3i  i  .pbi  31 1  ,aiilb(  31)  ,gi(  ii )  ,goi  3i ) 

V0«  1.138 

PO-  .101325 

A-  .1045 

B-  984.4-2.444*T 

IP  (P.LE.480)  go  to  50 

PH-. 101325«P 

Fx-  (a*pai/(a«P0) 

DO-LOG (PX I 
VLBZ-  V0/(1«A*ODI 
VSBZ-  .84a4*V0 
RETUBM 

SO  PBIITT*,  'EBBOB  *l»0  SOLID  EXISTS” 

BETURM 


c 

e 

SUBROUTINE  ESC4P 
c 

c  This  subroutlna  calculatas  tha  voluaa  of  cyclohaxana 

e  givan  a  praasura  and  taaparatura.  It  uaaa  publishod 


r 

r 

s 

I 


i 

I 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


'•V-' ' 


data  ftoa  Nohstn  Tasllai  a  doctotate  disettacion  1969, 
Tha  aquation  uaad  la: 

(  vo-v)/  Vo-  C  loqlO  (  B*P  )/  B 
(PI-  kp/caJ 
(81-  kp.'cal 
p.  41 


COmON/C/  P ,  VCPO ,  VCP ,  VOBZ ,  VLBZ  ,  VSBZ ,  V0C6 ,  VLC6 ,  TSH ,  MSB ,  POI ,  PBN 

CONHQH/0/  TH ,  KB  .NOKKt.B ,  CBCOST  .KOMKCY .  KNZSOt. .  HN2SPC.  CO ,  CS .  CIIL 

COmON/B/  T,CyCTH.S,ROU,e,WCOST,K 

COHMON/B/  PB<  31), PB'11),AMLR(11),CI(}1), 00(311 

B-  3002  -  6.a7«T 

C-  .2272 

VO-  1.284 

PM-  1.03127»p 

A-(aiPH)ya 

PRINT-. A 

VLC6-  VO*(1-C*LOG10(A) ) 

RETURN 

END 


SUBROUTINE  CRYSZ 

This  subcoutina  datcacainai  ainiaua  thall  and 
haad  thicknaat  of  tha  ptastuta  vatcel.  The  aquationa 
uaad  ara  ASMS  atandaeda  and  are  taken  fcoa  Max  S. 
pataca  and  Klaua  D.  Tianathaua,  Plant 
Oasiqn  and  Econoaioa  for  Chaaieal  Enoinaars, 

3rd.  Ed.,  p. 570-571. 

for  cylindrical  shall; 

t  -  r  (Z*».5  -1) 

Z  -  (SE  *  P)/'SE  -  Pi 

For  ellipaoidal  haada; 

t  -  2  Pr/(2SE  -  .2P) 

t  -  ainiaua  vail  thicxnaaa  (in) 

S  •  aaxiaua  allovable  workinq  stress  (ata) 

C  -  affiancy  of  joi''.ta  >a  a  fraction 
P  -  aaxiaua  intarnal  pressure  (psi) 

C  •  corrosion  allowance  (in) 


COMMON/  C/  P ,  VCPO .  VCP ,  VOBZ  ,  VLBZ  .  VSBZ .  V0C6 ,  VLC6 ,  TSH .  WSH ,  POI .  PRN 
COMHON/O/  TR  ,  WH  ,  MORE  LB .  CRCOST ,  MORKCY  .NNZSOL .  NN2SPC ,  CO  ,  C5 ,  CML 
COMMON  E/  T , C YCTM , S , ROU , E , NCOST , R 
COMMON/B/  rR( 31 ) , PRl 31 1 , ANLR( 31 ) .Gl (  31 ) .GO(  31 ) 

PI-  3.11159 

z-( ( s*et*pi/( (S»c)-P) 

TSH-  B«( (Z»».5)-l) 

TH-P»R»2/(2»S*E-.2»P) 

VC  -  VCPC-(4/3)*.0l619*PI*lR»*l) 

ALSH-  .3937»)VC*l55)/iPI*lB**2)) 

PRINT*, 'Z  L  ',Z, ALSH, VC 


VSH«  UTSB*I»)*»2-»**2)«»I*ALSH 

va-  («/31»H»(  tTH*«/2)*tTH*«)*»2-.S*««*3) 

HSR-  VSH*ROU*16.3S7/454 

HH-  VR*R00*16.3S7/454 

HT-  HSHtHH 

CRCOST-  WT«HCOST 

RSTUMI 

END 


stnaouTiNC  home 

In  this  aubroutin*  tha  cavttiabla  vork  is  csXculatad 
(ec  aach  at  tha  systaa  coaponantt.  It  is  assuaad  that  tha 
work  is  indapandant  of  coaposition.  Tha  aquations  of 
stats  usad  abova  for  datacaininq  tha  spacltlc  voluaa 
of  tha  liquids  ara  aqain  usad.  Tha  cavarsabla  work, 
tha  '.ntarqcal  of  Pdv,  is  ealculatad  by  tha  Ttapizoldal 
ruia  ovac  savacal  staps,  MSTEP. 

CORNON/C/  P . VCPO , VCP , VOSZ , VLBZ . VS»Z , V0C6 , VLC$ , TSH , HSH , POt , PRN 
CORRCN/O/  TH .MR , MORELS , CRCOST , WORKCT , WN2SOL , HN2SPC , CO . CS , CHL 
CORROS.'E/  T .  C YCTH ,  S ,  ROO .  B ,  WCOST,  R 
COHRON/B.'  PR(  31  ),PRi31l,AHLB(3i;,CI(31t,GO<  31) 

A-. 18«S 
C..2i7J 
NSTEP-30 
VOBZ-  1.138 
V0C6-  1.284 

07C6-  (V0CS-VLC6)/NSTCP 

BBZ-  9S4.6-2.44«*T 

8C«-  3002  -  4.e7*T 

PQ-.10132S 

00  lO.I-O.NSTEP 

VC«-  V0C«  -  I*0VC6 

CC«  IVOC«-VC6)/(VOC8*C) 

PC«MBH«  BC6»((10»»CC1-1)«1. 03327 
PC^.VVt.  (PC8NEH  ♦  PC«OtO)/2 
MORKZ-  NORR3  ♦  PC6AVE*OVC6 
PC60LD  •  PC6MBH 
CONTINUE 

po-i.q 

CPR-0 

OPS-0 

Opw  (P-P0)/NSTtP 

POLO-PO 

VO-VOBZ 

VSO-.9751 

DO  11,  J-l.NSTEP 

PNBH-  PO-  J*OR 

PAVE-  I  PNEN*PObO)/2 

IP' PAVE. LE. 1000)  GOTO  12 

cn-  (PAVE-POII  /<-2430) 

PRN-  .5»( (P«(1)*78*PR( J)«84»»(l-CN)-rR(ll«*4l/((l-CN)-CS) 

OPR-  PRN-PRO 

VSN-  V0BZ-t>PAVS*lS320)/101900) 

DVS  •  VSC-VSH 
IPiDVS.LE.0)DVS-0 

AA-  '.BBZ«PNEM)/1BBZ*P0) 
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VN«  VOBZ/(  l.*A*LOG(  AA)  ) 

OV-  -VM^VO 

NORRL>  OV«rAVt»(  rt(  1  I•78-PRII1 

IfOailT*  DPR ‘PAVE*.  0962 

VIORRS>  PRN*PAVt*OVS 

NORRl-  H0RR1><I0RRL*M0RRT*H0RRS 

VO-VM 

PRO-PRH 

VSO-VSN 

POLD-PNEN 

PRiirr*.'P  ave'.'horr  i*. 'product  rate* 

PRINT* .  PAVE ,  MORRl ,  PRN 
CONTINUE 

HORR4  •  HORRl  *MOa«2*E4*PR(  3) 

NORR4  •  NORE4*.4S4*. 0961118 

WORRCT-  HOBR4 

HORRLBa  H0RRCY/(PRN*2.2! 

RETURN 

END 


SUBROUTINE  N2REO 

This  subroutine  calculates  the  nitrogen  required 
by  the  crystalliier  for  the  bensene-n-cyclohexane  system. 
The  equation  of  state  used  is  the  real  qes  lav  vith 
the  coapressibillty  factors  found  by  estiaatinq  data 
to  a  straiqht  line  between  200  and  1000  ata.  The  data 
was  found  in  N.  Hooce,  Physical  Cheeistry,  4th  Ed. (1972), 
p.  20.  The  equqtion  is;  Z  »  1.3525C-3  P  ~  .75100.  There 
is  also  the  problea  of  the  nitraqen  disolvinq  in  the 
solution.  This  quanity  is  estinated  by  a  know  nele 
fraction  of  nitrcpen  in  benzene  at  25  oC  and  120  ata. 


COHHON/  C/  P ,  VCPO ,  VCP .  VOBZ  ,  VtBZ ,  VSBZ ,  V0C6 .  VtC6 ,  TSR,  NSH ,  PO I ,  PRN 
COHHON/D/  TH ,  NB ,  NORRLB .  CRCOST ,  NORRCV .  NN2SOI. ,  HN2SPC ,  70  ,  CS ,  CHL 
CONNON/E/  T ,  CrCTN ,  S  ,  ROi;,  E ,  NCOST,  R 
CONHON.'B/  PR(  311  .PR(  31)  ,AHLR<  31 )  ,Gt(  31 1  .CO(  31 ) 

SCI.  •  .02050 
VSPC  -  VCPO  -  VCP 

WN23PC  -  61.6*P*V8PC/((1.1525E-3*P*.38)*82.05*T) 

ODD-  PR<ll*78*rR(3)*64-PRN 

WN2SOI,  -  (DDO/(l-SOt)  )*SOL 

WN2  •  (NN2SOL  *  NN2SPC)/CVCTN 

aO(5)-WN2/26.02 

SI : 5)<NN2/28.02 

RETURN 

END 
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